Abstract. This is the first paper in a series of three where 
we present the first comprehensive inventory of sohd state 
emission bands observed in a sample of 17 oxygen-rich cir- 
cumstellar dust shells surrounding evolved stars. The data 
were taken with the Short and Long Wavelength Spectro- 
graphs on board of the Infrared Space Observatory (ISO) 
and cover the 2.4 to 195 /im wavelength range. The spec- 
tra show the presence of broad 10 and 18 fim bands that 
can be attributed to amorphous silicates. In addition, at 
least 49 narrow bands are found whose position and width 
. indicate they can be attributed to crystalline silicates. Al- 
' most all of these bands were not known before ISO. The 
^1;;^ ■ incredible richness of the crystalline silicate spectra ob- 

■ served by ISO allows detailed studies of the mineralogy 
of these dust shells, and is a telltale about the origin and 
evolution of the dust. We have measured the peak posi- 

, tions, widths and strengths of the individual, continuum 
subtracted bands. Based on these measurements, we were 
able to order the spectra in sequence of decreasing crys- 
talline silicate band strength. We found that the strength 
of the emission bands correlates with the geometry of the 
^— ^ ' circumstellar shell, as derived from direct imaging or in- 
, ferred from the shape of the spectral energy distribution. 
This naturally divides the sample into objects that show a 
disk- like geometry (strong crystalline silicate bands), and 
objects whose dust shell is characteristic of an outflow 
-«.^ , (weak crystalline silicate bands). All stars with the 33.6 

■ /im forsterite band stronger than 20 percent over contin- 
uum are disk sources. We define spectral regions (called 



> 



Oh- 



O ' complexes) where a concentration of emission bands is ev- 
^ ; ident, at 10, 18, 23, 28, 33, 40 and 60 /im. We derive 
^ ■ average shapes for these complexes and compare these to 
• • , the individual band shapes of the programme stars. In an 
. ^ ' Appendix, we provide detailed comments on the measured 
, band positions and strengths of individual sources. 
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1. Introduction 

Red giants and supergiants are characterized by low sur- 
face temperatures resulting in the presence of many differ- 
ent molecules in their atmosphere. These objects are also 
known to have dense stellar winds, presumably driven by 
a combination of pulsations and radiation pressure on the 
dust which forms in the cooling outflow. Since dust efh- 
ciently absorbes radiation at short wavelengths, the cen- 
tral star can easily become obscured and most of the lu- 
minosity of the star is re-radiated at mid-IR wavelengths. 
At these wavelengths, the most important ro-vibrational 
bands of abundant molecules can be found, and indeed 
many solid state bands from various dust components 
have been found using infrared spectrographs. The In- 
frared Space Observatory (ISO, Kessler et al. 1996) has 
allowed for the first time a comprehensive inventory of 
solid state bands in astrophysical objects with uninter- 
rupted wavelength coverage from 2.4 to 200 fim and with 
a spectral resolution which is well suited for the detection 
of solid state bands. 

We have undertaken detailed studies of the dust emis- 
sion and absorption spectra of evolved stars, ranging from 
Asymptotic Giant Branch (AGB) stars to Planetary Neb- 
Send offprint requests to: F.J. Molster: fmol- 
ster@so.estec.esa.nl 
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ulae (PNe) and to (post) Red Supergiants (RSG). Pre- 
liminary results of these studies have been published else- 
where and can be summarized as follows: Many oxygen- 
rich evolved stars have a surprisingly rich spectrum of solid 
state emission bands between 10 and 100 /im, dominated 
by both amorphous and crystalline silicates (e.g. Waters et 
al. 1996; Molster et al. 1999a; 1999b). The crystalhne sil- 
icates were not known to be present in dust shells around 
evolved stars before ISO was launched, and allow for the 
first time a mineralogical analysis of the dust composition 
around these objects. Some objects show a very high abun- 
dance of crystalline silicates (e.g. Molster et al. 2001a), 
which seems to be related to the geometry of the circum- 
stellar dust shell (Molster et al 1999a). Surprisingly, stars 
which were believed to have a carbon-rich dust chemistry 
also showed the presence of crystalline silicate emission, 
pointing to a complex chemical composition of the cir- 
cumstellar environment, possibly due to rapid changes in 
the chemistry of the stellar photosphere (e.g. Waters et 
al. 1998; Cohen et al. 1999, Molster et al. 2001a). These 
observations demonstrate the use of the crystalline silicate 
bands for a better understanding of the evolution of late 
type giants and supergiants. 

This paper is the first in a series of three where we 
present a detailed and comprehensive overview of the solid 
state emission bands in oxygen-rich dust shells surround- 
ing evolved stars and related objects. The purpose of these 
papers is to quantify as best as possible the presence and 
characteristics of the numerous new emission bands that 
have been discovered using the ISO data. In paper II (Mol- 
ster et al. 2002a) of this series, we describe the average 
band profiles of seven " complexes" that can be recognized 
in the combined Short- (de Graauw et al. 1996) and Long- 
(Clegg et al. 1996) wavelength spectrometers (hereafter 
referred to as SWS and LWS respectively) that were on 
board of ISO. Based on the strength of the crystalline sil- 
icate bands, we divide in paper II the sample of 17 stars 
into two groups. This division is also one which separates 
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Table 1. The programme stars and logbook of the observations used for this study. The boldface revolution numbers 

are those used for the final spectrum. 



Object 


Type 


Revolution 


AOT 


Tint 


Comment 








(AOT nr) 


(sec) 




IRAS09425-6040 


C-rich AGB star with O-rich dust 


084 


SWSOl 


1816 








254 


SWSOl 


6538 




NGC 6537 


Planetary Nebula, hot central star 


470 


SWSOl 


1912 


mispointed 






470 


LWSOl 


1318 








703 


SWSOl 


3454 








703 


LWSOl 


2230 




NGC 6302 


Planetary Nebula, hot central star 


094 


SWSOl 


6528 








479 


SWS06 


8532 


2.4 - 7.0 and 12.0 - 27.5Atm 






479 


SWS06 


12165 


7.0 - 12.0 and 27.5 - 45.2/im 






289-678 


LWSOl 


13243 


combination of rev 289, 482, 












489, 503, 510, 671, and 678 


MWC 922 


Peculiar object 


153 


SWSOl 


1834 


combined with rev 703 






478 


LWSOl 


1316 








703 


SWSOl 


1912 


combined with rev 153 


AC Her 


Binary post-AGB star 


106 


SWSOl 


1834 








471 


LWSOl 


1910 








520 


SWSOl 


6538 




HD 45677 


B[e] stax, nature unclear 


711 


SWSOl 


6538 




89 Her 


Binary post-AGB star 


082 


SWSOl 


1044 








336 


LWSOl 


1860 








518 


SWSOl 


6538 




MWC 300 


Evolved star, B supergiant 


516 


SWSOl 


3454 




Vy 2-2 


Proto-planetary nebula 


320 


SWSOl 


1140 








547 


LWSOl 


1318 


« 7 arcscc mispointed 


HD 44179 


Red Rectangle; binary post-AGB star 


702 


SWS06 


1174 


31.4 - 35.lMm 






702 


SWS06 


856 


19.5 - 25/im 






702 


SWSOl 


6538 








709 


LWSOl 


3428 








870 


SWS06 


8406 


12 - 19.5^im 


HD 161796 


Post-AGB star 


071 


SWSOl 


1044 








080 


LWSOl 


1554 








342 


SWSOl 


1912 








521 


SWS06 


1744 


29.0 - 45.2nm 


OH 26.5+0.6 


OH/IR star, high mass loss rate 


330 


SWSOl 


1912 








330 


LWSOl 


1268 








340 


LWSOl 


822 




Roberts 22 


post-AGB star, A supergiant 


084 


SWSOl 


1044 


mispointed 






103 


LWSOl 


478 








254 


SWSOl 


3454 


mispointed 


HD 179821 


post-AGB or post-RSG star 


113 


SWSOl 


1834 








520 


SWSOl 


6538 








319 


LWSOl 


1266 




AFGL 4106 


post-RSG, binary 


060 


SWSOl 


1130 








104 


SWSOl 


1834 








104 


LWSOl 


476 








249 


SWSOl 


3454 




NML Cyg 


Red Supergiant, high mass loss rate 


052 


SWSOl 


6544 








342 


SWSOl 


1140 








530 


SWS06 


1688 


29.5 - 45.2//m 






555 


LWSOl 


2798 








741 


SWSOl 


1140 




IRC+ 10420 


post-RSG, A supergiant 


128 


SWSOl 


3462 








316 


SWS06 


1718 


29.3 - 44.7/im 






724 


LWSOl 


3430 
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objects with a highly flattened dust distribution (referred 
to here as " disk" sources) from those with a more spheri- 
cal distribution of dust (the non-disk or spherical outflow 
sources, hereafter referred to as "outflow" sources). In the 
present study, we present the 17 programme stars, we give 
an overview of the individual spectra, and we quantify the 
position and strength of the bands. In Paper III (Molster 
et al. 2002b) we will investigate several trends in the spec- 
trum and correlate them with other information available 
about these sources. 

This paper is organized as follows: Sect. 2 presents 
the sample, the observations and data reduction; Sect. 3 
describes the nature of the individual sources, and the 
shape of the complexes as compared to the mean. Sect 4 
summarizes the results of this study. In Appendix A we 
present the band strength data and some more detailed 
discussion about individual spectral features with respect 
to their reliability and blending. 

2. The observations and data reduction 

2.1. The sample 

In order to get as broad as possible an overview of the 
0-rich dust features around evolved stars, we aimed for a 
sample covering the evolution of stars from the Asymp- 
totic Giant Branch (AGB) to the Planetary Nebula (PN) 
phase, and also including some massive (post-) Red Su- 
pergiants (RSG). Several objects were also included whose 
evolutionary status is unclear (e.g. MWC 922, MWC 300 
and HD 45677). We inspected the spectra for the pres- 
ence of crystalline silicate bands. Since we are interested 
in (sometimes weak) bands on top of a continuum, we se- 
lected bright objects with fluxes higher than about 20 Jy 
in the 20 to 45 /xm region. We preferentially included ob- 
jects for which both an SWS and an LWS spectrum are 
available. We finally arrived at a sample of 17 objects, 
listed in Table |l]. 

In several cases more than one observation was avail- 
able. The boldface revolution numbers in Table |^ indicate 
which observations were used for the final spectrum. The 
other observations were less accurately reduced and only 
used as a reference in case of doubt in the main spectrum. 

2.2. Data reduction 

2.2.1. SWS and LWS data reduction 

For a description of the SWS instrument and its operating 
modes we refer to de Graauw et al. (1996). The spectra 
used in this study were reduced using version 7.0 of the 
SWS offline processing software. For a description of flux 
and wavelength calibration procedures, see Schaeidt et al. 
(1996) and Valentijn et al. (1996). We used standard pro- 
cedures for flatfielding, sigma-clipping and rebinning of 
the signals from the 12 detectors We removed the main 
fringes in SWS band 3 (12.0 - 29.5 ^lm). The final reso- 



lution (A/AA) of the spectrum depends on the observing 
mode and the quality of the spectrum, but was typically 
around 250 for a speed 1, 300 for a speed 2, 500 for a 
speed 3, 750 for a speed 4 AOTOl and 1500 for an AOT06 
observation. 

The LWS instrument and operating modes are de- 
scribed in Clegg et al. (1996). For a description of the flux 
and wavelength calibration procedures see Swinyard et al 
(1996). The LWS spectra were reduced with the LWS of- 
fline processing software version 7.0, and further processed 
using ISAP. The reduced and rebinned LWS spectra of 
NGC 6302 used in this study were taken from Molster 
et al. (2001b); those of OH 26.5-1-0.6 from Sylvester et al. 
(1999); those of AFGL 4106 were taken from Molster et al. 
(1999b); for HD 44179 from Hony et al. (in preparation), 
for HD 161796 from (Hoogzaad et al., in preparation), 
and for MWC 922 from Sylvester (private communica- 
tion). We verified that these LWS spectra were reduced 
in a homogeneous way. The other LWS spectra were only 
sigma-clipped and rebinned and mainly used as a reference 
for the adjustment of SWS band 4 and the placement of 
the continuum (see section ^.2.2| ). 

To match the individual sub-bands with each other, we 
used multiplications, when we expected the flux calibra- 
tion to be the most important source of error, and a linear 
shift when we expect the the dark curents to be the main 
source for the discrepancy between the different bands. In 
all cases we tried to minimize the necessary shifts. 

2.2.2. Determination of an underlying contimmm 

In order to enhance the visibility of the different features 
we have defined a continuum for the spectra in Fig. |^ 
and Fig. |[ We have used an eye-ball spline-fit continuum, 
maximizing the continuum and still be smooth (no sudden 
changes in the slope), both in F^, and F\. We emphasize 
that this continuum has not necessarily a physical mean- 
ing but is only used here to enhance the sharp features 
on top of the spectrum. In principle, the strength of the 
features could be underestimated and there is a possibility 
that very broad features are treated as continuum. When- 
ever possible we tried to use the whole wavelength range 
(SWS -I- LWS) to determine the placement of the contin- 
uum. We have purposely not used (modified) blackbody- 
fits, because most spectra could not satisfactory be fit- 
ted by one (modified) blackbody (BB). This would have 
forced us to use multiple (modified) blackbodies, which 
would make it as arbitrary as the eye-ball fit and we had 
the feeling that we could be more consistent between the 
different spectra bye the eye-ball method than with BB- 
fitting. Finally, the BB-fitting results in a quite artificial 
physical parameter, the temperature, which is not always 
directly related to the real temperature of the dust species. 
In order to prevent confusion about this parameter, we 
prefered the meaningless eye-ball continuum. 



Fig. 1. The spectra of the disk sources (solid line) and their eye-ball spline-fit coritiriuurn (dashed line). The jumps in the 
spectrum of Roberts 22 are due to a mispointirig of the satellite. The straight lines after 45/im in the spectra of IRAS09425- 
6040, HD 45677 and MWC 300 are the connections with the IRAS datapoints. 



In two sources, Roberts 22 and Vy 2-2, flux jumps were 

found due to a pointing offset of ISO. Pointing offsets re- 
sult in flux jumps due to the change in aperture size of 
SWS with wavelength (see De Graauw et al. (1996)). In 
these cases wc have fitted separate continua to each part 
of the spectrum, where we took into account the relative 
slope in the other parts of the spectrum. However, this 
may affect the shape of the complexes whose wavelength 
coverage extends over sub-bands with different SWS aper- 
tures. 



The Planetary Nebulae NGC 6537 and NGC 6302 have 
strong emission lines that hamper the deternnnation of the 
dust emission bands. We decided to remove most of these 
by subtracting a Gaussian flt to the (unresolved) lines. For 
very strong lines (some could be 100 times the continuum 
level) the subtraction procedure resulted in excessive noise 
after removal of the Gaussians. In these cases, we removed 
the noisy part of the spectrum. 



Fig. 2. The spectra of the outflow sources (sohd hne) and their eye-ball-spUne-fit continuum (dashed line). The jumps in the 
spectrum of Vy 2-2 are due to a mispointing of the satellite. 



2.2.3. The "final spectra": disk and outflow sources 

After some study of the stars in our sample, it became ap- 
parent that there is a wide range in band strengths of the 
different sources. We have accordingly ordered the spec- 
tra in sequence of decreasing strength of the 33.6 /im band 
with respect to its local continuum. As we will show later 
(see Sect. ||) this naturally divides the spectra into ob- 
jects that have a disk-like distribution of dust, and ob- 
jects that have a more spherical dust distribution. We will 
use this division into "disk" and "outflow" sources in the 
remainder of this paper, and will group our sources ac- 
cordingly. In Fig. |l| and Fig. || the final spectra are shown, 



and Figs. |3| and ^ show the continuum subtracted spectra, 
in both cases divided into disk and outflow sources. 

2.3. Definition of complexes 

A first inspection of the continuum subtracted spectra of 
our programme stars indicates that there is a multitude 
of emission bands. These bands are not evenly spaced in 
the spectrum, but tend to concentrate in a number of 
wavelength regions. Within these fairly narrow wavelength 
ranges, severe blending hampers the measurement of indi- 
vidual components. As an example we show in Fig. |^ the 



Fig. 4. The continuum subtracted spectra of the outflow 
sources. 

OH 26.5+0.6 and NML Cyg, because they show the 10 
/im amorphous silicate band in strong absorption. 

The source OH 26.5+0.6 has also been excluded from 
the average spectra of the 18 and 23 micron (outflow) com- 
plexes because the spectrum is in absorption. NML Cyg 
shows amorphous silicates in absorption in the 18 micron 
complex, but the crystalline silicates is already in emis- 
sion. To avoid problems we have also excluded this object 
from the average outlfow spectrum for the 18 micron com- 
plex. 

In Paper II an extensive discussion will be presented 
on the properties of the mean complexes and bands, for 
both outflow and for disk sources. Here we will compare 
the individual complex or band positions and shapes with 
the mean. In Paper II we will also present average proper- 
ties and identify the bands using laboratory data of sev- 
eral materials. Most of the bands can be identified with 
forsterite (Mg2Si04; the Mg-rich endmember of the crys- 
talline olivine series) and enstatite (MgSiOa; the Mg-rich 
end member of the crystalline pyroxene series). About 20 
percent of the bands still lack an identification. 

2.5. Measurement of solid state features 

In order to get a consistent view of the strength and po- 
sition of the features found in our spectra, we have fitted 
Gaussians to features, or multiple Gaussians in the case 
of complexes, with the ISAP-line fit programme. The fits 



Fig. 5. The continuum subtracted spectrum of HD 179821. We indicated the range of the different complexes. 



Table 2. The definition of the different complexes 



complex 


^min 


Amax 


name 


/im 


Hm 


10 


7 


13 


18 


15 


20 


23 


22 


25.5 


28 


26.5 


31.5 


33 


31.5 


37.5 


40 


38 


45.5 


60 


50 


72 



have been applied to the original, i.e. non-continuum sub- 
tracted, spectriim in Fx. In this way wc prevent very broad 
features, extending over more than one complex, from con- 
tributing to the individual bands. We determined a local 
continuum, often a third order polynomial, for each spec- 
tral section and estimated the error on the measured wave- 
length, FWHM, peak to continuum ratio and integrated 
band flux for these features by multiple independent fits to 
the (local continuum subtracted) spectrum. For this pur- 



pose we varied the defined continuum, since this is likely 
to have the largest influence on the derived strength of 
the features. The drawback of this method is that the un- 
certainties in the spectrum are not taken into account. 
Therefore in a few low S/N cases the derived 1 a errors 
for the wavelength position (see Appendix A), might be 
too low and should be a factor 2 to 3 higher. Still, for most 
cases the wavelength determination is indeed quite accu- 
rate. Sometimes the error in the strength of the features 
may have been overestimated, resulting in band strengths 
that arc less than 3 sigma over the noise. However, care- 
ful inspection of the individual cases convinced us of their 
reality. Finally, we note that non-Gaussian shapes occur 
for certain bands, causing a systematic but reproducible 
error. 

3. Description of the programme stars 

Before discussing the bands and complexes, it is useful to 

summarize the properties of our programme stars. This 
is relevant when one tries to correlate the observed dust 
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Table 3. Overview of the different features soon in our sample (l=IRAS09425-6040; 2=NGC 6537; 3=NGC 6302; 4=MWC 922; 
5=AC Her; 6=HD 45677; 7=89 Her; 8=MWC 300; 9=Vy 2-2; 10=HD 44179; 11=HD 161796; 12= OH 26.5+0.6; 13=Roberts 22; 
14=HD 179821; 15=AFGL 4106; 16=NML Cyg; 17=IRC+ 10420). The wavelengths are in /xm. A '?' indicates an uncertain 
detection. Blends are indicated by '}'. ^ indicates that this features is only found in 1 star. 17.5 and 17.5b indicate respectively 
a narrow (enstatite) band and a broad (amorphous silicate) band. The features at 11.05, 13.5 and 32.97 nm are instrumental 
artifacts. The status of the 13.8 and 14.2 micron features is unclear, they also seem to suffer from instrumental artifacts, but a 
contribution from onstatito cannot be excluded. 



8.3 8.20 

9.14 \ 9.12 9.17 9.15 9.12 

9.45 J-9.3a 945 945 945 945 

9.8 9.84a 9.83 

10.1 9.97 9.59 10.0 9.9a 10.05 10.4 9.9a 10.57 10.58 9.94a 10.6 

10.7 10.6? 10.60a 10.90 10.57 10.68 10.8? 10.7? 
11.05 11.06 11.06 11.06 11.07 11.06 11.04a 

11.4 11.37 
11.9' 11.90 

13.5 13.60 13.55 13.54 13.54 13.51 13.53 13.58 13.4? 13.40 \ 

13.8 13.90 13.84 13.77 13.87 13.74 J13.6 

14.2 14.2? 14.22 14.19 14.3 14.20 14.15 14.19 14.16 14.20 14.24 14.17 
15.2 15.13 15.24 15.35 15.37 15.0? 15.02 15.42a 

15.9 15.79 15.92 15.94 15.69 15.87 15.88 16.06a 
16.2 16.19 16.17 16.10 16.19 16.18 16.37a 



16.50 



18.09 


18 


11 


19.11 


18 


97 


19.60 


19 


56 


20.66 


20 


67 


21.6 


21 


62 




22 


38 


22.98 


22 


97 


23.63 


23 


65 


24.50 


24 


54 


25.12 


25 


00 


26.29 






26.96 








27 


54 


^28.0 


28 


04 


28.70 






29.62 


29 


65 


30. .55 


30 


51 




31 


12 


32. .34 


32 


25 


32.83 


32 


78 


33.70 


33 


71 


34.36 


34 


24 




34 


99 


j-35.1 


35 


82 


36.56 


36 


56 


39.96 


39 


9 


40.53 


40 


48 


41.58 


41 


55 


42.99 


42 


74 


43.74 


43 


55 


44.67 


44 


58 


47.70 


47 


62 


48.64 


48 


49 


51.40 


52 


1 


j-62.4 


}-63.5 


68.79 


68 


96 


90.7 


91 


1 



8.31 


8.40 






9.17 


9.15 






9.45 


9.46 






9.81 


9.8? 






9.59 


10.0 


9.9 


a 


10.57 


10.68 






11.06 


11.06 






11.50 


11.33 






13.54 


13.51 


13 


53 


13.77 


13.87 


13 


74 


14.3 


14.20 


14 


15 


15.37 


15.0? 






15.94 


15.69 








16.10 






j-17.0 


}l6.9 


16 


98 


17.57 




17 


61 


18.5 


18.0 


17 


Oa 


18.06 


17.95 


17 


96 


18.88 


18.49 


18 


43 


19.43 


19.36 


19 


36 


20.67 




20 


66 



}l6.1 j-16.1 j-16.1 

|-16.9 }l7.0 }l7.0 J70I }l6.8 }l7.0 |-16.9 iggg }l7.1 }l7.1a }l6.9 }l 

17. 5n 17.43 17.54 17.57 17.61 17.51a 

17.5b 16.8a ?? 18.5 18.0 17.0a ?? 17.4 17.78a 18.0 17.7 17.70 

1.0 18.09 18.09 18.11 18.15 17.90 18.06 17.95 17.96 18.03 18.04 18.16a \ 18.01 

18.9 19.17 19.11 18.97 18.94 18.75 18.88 18.49 18.43 18.97 18.84 18.99a j-18.3 igf 

19.5 19.63 19.60 19.56 19.39 19.64 19.43 19.36 19.36 19.65 19.49 19.53a 19.75 19.5 

20.7 20.69 20.66 20.67 20.64 20.84 20.67 20.66 20.58 20.77 20.54 20.57a 20.70 20.65 20.63 

21.5 21.49 21.6 21.62 21.65 21.42 21.59 21.35 21.59 21.52 21.42 21.43 

22.4 22.43 22.38 22.37 22.41 1 22.51a 22.49 22.26 
22.69' 22 
23.0 23.03 22.98 22.97 22.88 1 

23.7 23.65 23.63 23.65 23.45 J23.5 23.67 23.67 23.72 IjS 3 ^3.81 23.65 J-23.30a 
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spectrum with other properties of the object. We have 
searched for correlations between the relative (to the lo- 
cal continuum) strength of the 33.6 /xm band (identified 
with forsterite, Waters et al. 1996), which is detected in 
all programme stars, and other properties of the star or 
its circumstellar shell. No obvious correlation could e.g. 
be found with the temperature of the star, or with the 
colour temperature of the underlying continuum (we will 
return to these points in Paper III). However, a correlation 
between the geometry of the circumstellar shell and the 
(relative) strength of the 33.6 nm feature exists (Molster 



et al. 1999a). Objects with strong crystalline silicate emis- 
sion bands have a highly flattened circumstellar geometry. 
This geometry is either derived from direct imaging at op- 
tical or near-IR wavelengths (e.g. in the case of the Red 
Rectangle: Osterbart et al. 1997; Monnier et al. 1997), or 
it is inferred from the shape of the spectral energy distri- 
bution (SED; e.g. in the case of IRAS09425-6040, Molster 
et al. 2001a). 

Because of these differences we have divided the stars 

into disk and outflow sources and ordered them in decreas- 
ing 33.6 micron band strength. When discussing individual 
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stars below, present evidence for cither a disk or a more 
spherical geometry for the dust envelope. 

3.1. Disk sources 

3.1.1. IRAS09425-6040 

This carbon-rich AGB star is one of the most intriguing 
stars in our sample. The source is discussed extensively 
by Molster et al. (2001a). It was classified as a post- AGB 
object based on its IRAS colours. The ISO-SWS spectrum 
however revealed a J-type carbon-star at wavelengths be- 
low 15 /J,m (Molster et al. 2001a), while at longer wave- 
lengths highly crystalline 0-rich dust dominates the spec- 
trum. It has the highest contrast crystalline circumstellar 
dust observed so far in the complete ISO sample (Molster 
et al. 1999a). The oxygen-rich dust probably originates 
from a previous mass loss phase before the star turned 
C-rich. Somehow this 0-rich material was stored around 
the star. The IRAS 60 /zm measurement suggests that the 
spectral energy distribution is flat up to that wavelength. 
Both the shape of the energy distribution and the com- 
plex chemistry point to the presence of a disk instead of 
an ordinary outflow. Molster et al. (2001a) propose that 
IRAS09425-6040 is a binary and an evolutionary progen- 
itor of the Red Rectangle. 

The SWS spectrum of this source has already been 
shown by Molster et al (1999a;2001a). 

3.1.2. NGC 6537 

NGC 6537 is an extreme Type I planetary nebula (PN), 
which bears a lot of similarities to NGC 6302. It has a 
bipolar outflow structure, which is probably caused by a 
disk (Cuesta ct al. 1995). He and N are overabundant 
while C is severely underabundant, the typical character- 
istic of Type I PN, but for both NGC 6537 and NGC 6302 
these abundances are more extreme than in other Type I 
PN. Both NGC 6537 and NGC 6302 have very high ex- 
citation hues, e.g. [Si VI] (Ashley & Hyland, 1988) which 
are not foimd in other PNc. NGC 6537 is probably the 
evolutionary product of the most massive stars which end 
as a white dwarf and may have had an initial main se- 
quence mass in the range 5-7 Mq. Shocks arc expected 
to be present due to the interaction of the fast wind and 
the disk. The infrared spectrum shows evidence for both 
carbon-rich dust (PAHs) and oxygen-rich dust (silicates). 

3.1.3. NGC 6302 

NGC 6302 or the 'Butterfly Nebula' is one of the infrared 

brightest PN. It has one of the hottest central stars known 
with Teff = 250000 K (Casassus et al. 2000). Together 
with the high abundance of nitrogen and helium in the 
nebula (Aller et al. 1981) this points to a massive progen- 
itor. It is a Type I PN, with a highly bipolar outflow and a 



thick dusty 'disk'. The spectrum shows highly excited for- 
bidden emission lines (up to [Si VII] for NGC 6302). Roche 
& Aitken (1986) already detected the PAH-features at 8.6 
and 11.3 micron, and Cohen et al. (1989) presented evi- 
dence for the 6.2 and 7.7 micron PAH feature. Both obser- 
vations indicate the presence of carbon-rich material. On 
the other hand Payne et al. (1988) detected an OH-maser, 
normally associated with oxygen-rich environments. This 
dichotomy in the dust is also present in our ISO spec- 
trum, where we find the PAH-features and the crystalline 
silicates. The presence of an equatorial disk-like structure 
has long been known (Meaburn & Walsh 1980; Lester & 
Dinerstein 1984; Rodriguez et al 1985) and can also be 
seen as a dark (dust) lane on images of this PN. The 30 
to 45 /xm part of the spectrum has first been published by 
Waters et al. (1996), Beintema (1998) showed the com- 
plete SWS spectrum while the LWS spectrum was shown 
by Barlow (1998). 

The fufl ISO spectrum of NGC 6302 has been discussed 
by Molster et al. (2001b). 

3.1.4. MWC 922 

The evolutionary status of MWC 922 is unclear, it was 
placed among the unclassified B[e] stars by Lamers ct al. 
(1998). Both a pre-main sequence (The et al 1994) as well 
as a post main sequence status (Voors 1999) have been 
suggested. The distance to this object is unknown. The 
star is projected on the Ser OBI association, which is 
at 1.7 kpc. However, there is no independent evidence 
that MWC 922 should be associated with Ser OBI. Si- 
mon & Dyck (1977) found an infrared excess at 20 and 
25 /um. Meixner et al (1999) observed this object at 8.2 and 
12.2 /xm and marginally resolved it. Its broad spectral en- 
ergy distribution and rather high mm continuum flux are 
not compatible with a spherically symmetric (continuous) 
dusty outflow, and we classify the star as disk-like. Spectra 
taken with the Kuiper Airborne Observatory (KAO) data 
show strong PAH emission (Cohen et al. 1989), which usu- 
ally is found in carbon-rich environments. The ISO spec- 
trum shows that oxygen-rich material is also present. 

The continuum divided LWS part of the spectrum was 
already shown by Barlow (1998), while the SWS part was 
shown by Voors (1999). 

3.1.5. AC Her 

AC Herculis is an RV Tauri star, with a stable pulsation 
period of 75.47 days (Zsoldos 1993). The star was foimd 
to be a binary with a period of 1194 days (van Winckel 
et al. 1998). The narrow velocity width of the CO rota- 
tional line emission (Jura et al. 1995) suggests Keplerian 
rotation in a disk rather than outflow. Also, very strong 
mm continuum flux (van der Veen et al. 1994) indicates 
the presence of large grains, which suggests a long stor- 
age time and therefore a (planet forming?) dust disk. The 
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infrared spectrum of the binary star AC Her is strikingly 
similar to the spectrum of the comet Hale Bopp (Molster 
et al. 1999a). The SWS spectrum longwards of 7/um has 
already been shown by Molster et al (1999a) and is briefly 
discussed by van Winckel et al. (1998). Recently, the evi- 
dence for a disk has been provided by Jura et al. (2000). 
They found a dust ring with a radius of approximately 
300 AU in an image taken at a wavelength of 18.7 /xm. 

3.1.6. HD 45677 

HD 45677 is a well studied B2 star whose evolutionary sta- 
tus is still unclear. A pre- main sequence nature has often 
been suggested for this star, however its isolated position 
(HD 45677 is not associated with nebulosity) and the ab- 
sence of Algol-like variations and 'blueing' effect make this 
questionable. Lamers et al (1998) also discuss this star and 
place it amongst the 'unclassified' B[e] stars. They propose 
that it is an extreme example of a classical Be star. Po- 
larization measurements show that the circumstellar ma- 
terial is located in a disk (Schulte-Ladbeck et al. 1992). 
Also, the strong absorption cores of the Na I D and Ca 
II K lines indicate the presence of an optically thick disk 
at zero velocity (de Winter & van den Ancker 1997). This 
disk was already present before 1950 when a significant 
disruption took place. After this event, the star dimmed 
up to » 2 mag in 1981 without significant change in the 
colours. This is most easily explained by the production 
of large (= grey) particles (de Winter & van den Ancker 
1997). 

The spectrum of this star has already been published 
by Voors (1999) and Malfait (1999), who also modeled it. 

3.1.7. 89 Her 

89 Her is a high galactic latitude F2 Ibe supergiant with a 
temperature of about 6500 K (Waters et al. 1993). Waters 
et al. (1993) confirmed the binary hypothesis of Arrelano 
Ferro (1984), finding an orbital period of 288.4 days. The 
CO(l-O) and C0(2-l) line observations show a narrow 
(< 1 km/sec) central spike on top of a broader (w 8 
km/sec) weak component (Likkel et al. 1991). This pro- 
file is very difii'erent from the profiles seen from detached 
shells, and may represent gas orbiting in a fiattened disk- 
like structure. Several other arguments for the presence 
of a disk, such as the lack of energy balance between the 
UV and IR, are given by Waters et al. (1993). Alcolea & 
Bujarrabal (1995) imaged 89 Her in CO(l-O) and found 
an outer shell, which likely originates from a heavy mass 
loss period experienced by the star in the past. This erup- 
tion of mass was probably triggered by the companion 
by a Roche lobe overflow or even by a common envelope 
phase, and ended the AGB evolution of 89 Her. It is likely 
that during this period also the disk was formed. 



3.1.8. MWC 300 

MWC 300 is classified as a Bl Ia+ hypergiant by Wolf & 
Stahl (1985). With a luminosity of « 5 x lO^Lg the 
star is at 15.5 kpc and about 500 pc above the Galactic 
plane. Henning et al. (1994) detected MWC 300 at sub- 
mm wavelengths and pointed out that for the derived dis- 
tance of 15.5 kpc (Wolf & Stahl, 1985) the total dust mass 
in the circumstellar envelope would be in the order of 300 
Mq, which seems unrealistically high. Different attempts 
were made to spatially resolve this object (Skinner et al. 
1993; Ageorges et al. 1997; Leinert et al. 1997; Pirzkal et al 
1997), but they were all unsuccessful. This indicates that 
the dust is indeed circumstellar and not associated with 
nearby nebulosity. A more realistic circumstellar dust and 
gas mass of a few would result in a much smaller dis- 
tance (an order of magnitude) and therefore a luminosity 
of about 10" LQ.This new luminosity together with its ex- 
tended atmosphere characteristics would classify this star 
as an evolved low mass star (Voors 1999). Winckler & Wolf 
(1989) and Hamann and Persson (1989) both argue that 
MWC 300 is surrounded by a (slowly expanding) disk-like 
structure and a low density wind near the poles. 

3.1.9. HD 44179 

HD 44179 is the central star in the X-shaped Red Rectan- 
gle nebula (Cohen et al. 1975). It is an A type supergiant 
with an effective temperature of about 7500 K in a binary 
system (P=318 days), surrounded by a circumbinary disk 
(Waelkens et al. 1996), which we see (almost) edge on. The 
central star is heavily obscured by this disk and only seen 
in reflection by scattering lobes below and above the plane 
of the disk. The optically thick disk has been imaged with 
high resolution in the optical and near-IR by e.g. Roddier 
et al. (1995) and Osterbart et al. (1997). The CO(l-O) 
and C0(2-l) show very narrow (« 5 km/sec) line emission 
(Jura et al. 1995), atypical for detached AGB remnants. 
The mm and cm continuum flux of this source is rather 
high and suggests the presence of large (mm-sized) grains 
(Jura et al 1997). These grains are likely to be formed 
in the long-lived circum-binary dust disk. The similarities 
with the disks around young stars lead to speculations 
about possible planet formation around this evolved star 
(Waters et al. 1998). The detection of a mysterious dust 
clump around HD 44179 by Jura & Turner (1998) feeds 
this interesting speculation. The star has both carbon- 
rich, as evidenced by the PAH-features, and oxygen-rich 
dust, as evidenced by the crystalline silicate features. The 
PAHs are predominantly present in the scattering lobes, 
while the crystalline silicates are expected to be present 
in the disk (Waters et al. 1998). 

The SWS part of the spectrum was already shown by 
Waters et al. (1998). 
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3.1.10. Roberts 22 

Roberts 22 is a bipolar reflection nebula, whose evolution- 
ary status has recently been determined as being post- 
AGB (see e.g. Sahai et al. 1999). Allen et al. (1980) found 
that the central star is completely obscured by a central 
dust lane and they determined the spectral type of this 
object from the two reflection lobes, which gave identical 
spectra (A2 I). At a distance of 2 kpc (Allen et al. 1980; 
Sahai et al. 1999) its total luminosity is « 3 x 10** Lq. 
Prom IR and the CO emission line data the progenitor 
(AGB) mass- loss rate was about 10"'* M0/yr. Roberts 22 
has a time- variable OH maser (Allen et al. 1980), which is 
mainly located in the central waist, but is also seen in the 
northern and southern scattering lobes (Sahai et al. 1999). 
The velocity distribution of the OH masers might be in- 
terpreted as a rotating disk, which is seen almost edge-on. 
The fact that the OH is also seen in the the lobes suggests 
that the disk is being disrupted, probably due to the fast 
wind (450 km/sec) seen in Ha (also time variable) arising 
from the central star. The spectrum of Roberts 22 also 
shows both PAH-features as well as the infrared features 
of crystalline silicates. 

Part of the spectrum was shown by Molster et al. 
(1997). 

3.2. Outflow sources 

3.2.1. Vy 2-2 

Vy 2-2 is classifled as a compact planetary nebula sur- 
rounded by a fossil molecular envelope from the progeni- 
tor AGB star (Jewell et al. 1985 and references therein). 
Lamers et al. (1998) classified this star as a compact plan- 
etary nebula B[e] star. This very young PN has both an 
ionized zone and a neutral, molecular cloud. The ionized 
nebula has been resolved as a thin shell extending to ~ 0.5 
both at 15 Ghz (Seaquist & Davis 1983) and in Ha (Sa- 
hai & Trauger 1998). The inner radius was estimated to 
be 0.2. Only the blue shifted OH maser component is de- 
tected in Vy 2-2. This is not in contrast with a homo- 
geneous spherically symmetric outflow, because that sce- 
nario predicts that the red shifted component is absorbed 
(Seaquist & Davis 1983). Although the quality of the spec- 
trum is low, due to mispointing, we decided to keep it in 
our sample since it was the only young 0-rich PN. 

3.2.2. HD 161796 

HD 161796 is a high galactic latitude F3 lb supergiant. 
Skinner et al (1994) determined a distance of 1.2 kpc and 
therefore a luminosity of 3600 Lq. Its photosphere has 
low metallicity and an enhanced nitrogen abundance, im- 
plying that it is an evolved Population H object (Luck 
et al. 1990). The shape of the very strong CO emission 
found by Likkel et al. (1991) resembles the profile of mass- 
losing AGB stars and confirms the post-AGB nature of 



this object. The expansion velocity determined from this 
line is about 12 km/sec, a typical value for AGB stars. 
Skinner et al. (1994) resolved the envelope around this 
source in the mid-IR (10.5 and 12.5 /im). They foimd 
an expanding dusty equatorial toroid in a final phase of 
strongly enhanced, equatorially concentrated mass loss, 
which stopped about 240 years ago. During this mass- 
loss burst the mass- loss rate was about 3 x 10~^ M0/yr. 
Meixner et al. (1999) confirmed these results. The con- 
tinuum divided LWS spectrum was already presented by 
Barlow (1998). 

3.2.3. OH 26.5+0.6 

OH 26.5+0.6 is an extreme OH/IR star, which shows ev- 
idence of two mass-loss regimes: a superwind phase in 
which the mass-loss rate is ~ 10""' M0/yr which started 
about 200 years ago, and an earlier AGB phase with a 
mass-loss rate of about w 10~^ Mq/jt (Justtanont et al. 
1994, 1996a). The transition between these two phases 
was probably very short {At < 150 yr). The total mass 
lost during the superwind phase has been estimated to be 
« 0.1 Mq (Justtanont et al. 1996a). Infrared speckle in- 
terferometry at 9.7 /xm (in the silicate feature) gives an 
angular size for the circumstellar dust shell of 0.50 ±0.02, 
while the angular size of the dust shell outside this fea- 
ture (at 8.7 /Ltm) is less than 0'.2 (Fix & Cobb 1988). This 
difference is caused by the enhanced opacity in the 10 
lim silicate feature, therefore the dust seen at these wave- 
lengths is located at larger radial distances and is cooler 
than the dust seen on either side of the feature. 

The spectrum of OH 26.5+0.6 has already been shown 
by Sylvester et al. (1999). 

3.2.4. HD 179821 

HD 179821 is a G5 la supergiant at a distance of about 6 
kpc (Zuckerman & Dyck 1986; Hawkins et al. 1995; Jura 
& Werner 1999). At this distance the star would have a 
luminosity of about 3.1 x 10''' L0, far above the AGB lumi- 
nosity limit. If the star is indeed massive, the origin of the 
dust envelope was probably the Red Supergiant phase. In 
that case, the star may evolve to the Wolf-Rayet phase, 
before exploding as a supernova. The dusty envelope has 
been resolved at MIR wavelengths (Hawkins et al, 1995; 
Jura & Werner 1999). A ring- like structure was foimd with 
an inner radius of 1.75 (1.6 x 10*^ cm at a distance of 6 
kpc). This implies that the mass loss burst has stopped 
about 1500 years ago. In CO (Bujarrabal et al. 1992) and 
NIR scattered light (Kastner & Weintraub 1995) the dust 
envelope was found to extend to at least 18 , indicating 
that the mass loss burst lasted for at least 6000 yr. With 
a derived gas mass loss rate of the order of 10~^ M0/yr 
(Kastner & Weintraub 1995), about 6 Mq was lost by the 
star during this mass- loss burst. This value is similar to 
that found by Jura & Werner (1999) based on submillime- 
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tre measurements, if a gas-to-dust ratio of 200 is assumed. 
From infrared imaging (Jura & Werner 1999) and maps 
in CO (Bujarrabal et al. 1992), it is seen that the gas and 
dust distribution is not spherical. Part of the spectrum of 
this source was already shown by Waters et al. (1996). 

3.2.5. AFGL 4106 

AFGL 4106 has been discussed by Molster et al. (1999b), 
and is a high mass (15 to 20 Mq) binary with two almost 
equally luminous stars with temperatures of about 3750 K 
and 7250 K. The binary is located at a distance of about 
3.3 kpc. During about 4.3 x lO'^ years, the more massive 
star in the system (now the warmer of the two) had a gas 
mass-loss rate of about 9x10^'' M0 / yr, which gives a total 
expelled mass of about 4 Mq. This huge mass- loss burst 
ended about 450 years ago. Mid-IR imaging by Molster et 
al. (1999b) shows an asymmetric detached dust shell. 

3.2.6. NML Cyg 

NML Cyg is an M6 supergiant at a distance of 1.8 to 2 kpc 
(Morris & Jura 1983; Bowers et al. 1983) with a luminosity 
of « 5 X 10^ Lq, implying a main sequence mass of 50 Mq. 
The present-day mass loss rate is between 1.1 x 10~^ (from 
OH and IR; Netzer & Knapp 1987) and 1.8 x lO'^'Moyr"! 
(from CO; Knapp et al. 1982). The dust shell has been 
partially resolved by Fix & Cobb (1988^ in the N-band 
(10 fim). They found an outer radius of 0.37, which corre- 
sponds to about 10^^ cm. Monnier et al. (1997) found the 
same value for their outer dust shell, and detected a sec- 
ond dust shell, located inside the first one, at about 0.125. 
The outflow velocity determined from the 1612 MHz OH 
maser line is 27.7 km/sec (Bowers et al. 1983). In the H2O 
maser map an asymmetry is found at subarcsecond-scales 
(Richards et al. 1996). This is also seen in the OH maser 
maps on scales of a few arcsec, indicating that this asym- 
metry comes from the inside and is not due to external 
factors. Parts of the ISO spectrum were already shown by 
Justtanont et al (1996b) and Waters et al. (1996). 

3.2.7. IRC+10420 

IRC-l-10420 is an A5 Ia+ hypergiant with a large IR- 
excess. It is found to be continuously evolving from a 
cooler (F8 Ia+; Humphreys et al. 1973) to a hotter spec- 
tral type (A5 la"*"; Oudmaijer 1998 and references therein). 
IRC-l-10420 is at a distance of 5 ± 1 kpc which gives it a 
luminosity of about 5 x 10^ L0, just below the Humphreys- 
Davidson limit. From CO and OH measurements, an out- 
flow velocity of 40 km/sec has been derived (Nedoluha & 
Bowers 1992; Bachiller et al. 1988; Lewis et al. 1986). Kast- 
ner & Weintraub (1995) found from their near-IR polari- 
metric maps that the dust envelope extends to 9 , which 
implies a dynamical age of the dust in the outer layers 
of about 5000 years. The mass loss rate that created the 



dust shell was of the order of 10~^ M0/yr, and a total of 
5 M0 of gas and dust is present in the shell: a significant 
fraction of the 40 M0 which it probably had on the main 
sequence. Bowers (1984) observed this star in the 1612 
and 1667 MHz OH maser lines and found indications for 
multiple shells in OH extending out to 4 . This suggests 
that the mass loss went in bursts and was not constant 
in time. The circumstellar dust shell has been resolved in 
the mid-IR at 8.7, 9.8 and 20.6 /zm (Fix & Cobb 1988; 
Jones et al. 1993; Meixner et al. 1999). In these images an 
elliptical structure is found, which Oudmaijer et al. (1994) 
attribute to a bipolar outflow, which is beamed into our 
direction. 

4. Description of the complexes 

In this section we discuss the individual spectra with re- 
spect to the mean spectra and with respect to each other. 
In Fig. ^ to Fig. ^ we show the emission (and in a few 
cases the absorption) complexes of the individual stars and 
compare them with the mean spectra for every complex. 
In these figures we indicate the noise level and the wave- 
length (spread) of the individual features. The noise level 
is the mean 3a noise level, but might change with wave- 
length in one complex; e.g. due to a significant change 
in continuum level or to changes in SWS-band and there- 
fore detector characteristics. These effects are particularly 
present in the 10 and 28 micron complexes. Here, we limit 
the discussion to the featuress that we are confident are 
real. In Appendix A we discuss some more dubious detec- 
tions and artifacts. 

4-.1. The 10 micron complex 

The 10 micron complexes of the disk and outflow sources 
are shown in Fig s. I and §. 

4.1.1. Disk sources 

IRASO9425-6O4O: The 10 micron complex of IRAS09425- 
6040 is still dominated by the carbon star and its present 
outflow, therefore C-rich molecules are found in absorp- 
tion and prominent SiC emission is seen in the 10 micron 
complex. No clear evidence for crystalline material has 
been found. 

NGC 6537: The 10 micron complex of NGC 6537 is 
dominated by PAH features, and no evidence for (crys- 
talline) silicate features has been found. 

NGC 6302: The 10 micron complex of NGC 6302 is 
dominated by PAH features, and no evidence for (crys- 
talline) silicate features has been found. 

MWC 922: The 10 micron complex of MWC 922 is 
dominated by PAH features. Around 10 micron, absorp- 
tion features seem to be present, which line up very well 
with the position of crystalline silicate emission features 
in other stars. Since crystalline silicates are not abundant 



Fig. 6. The 10 micron complex for the disk sources, together with the mean disk spectrum (gray line). The thin gray vertical 
lines indicate the mean peak positions of the features found, while their range is indicated by the errorbar on top of the plot. 
The thick line in the upper left corner in each panel indicates the mean noise level. 



in the ISM, this must be circumstellar material, It is one 
of the few sources where crystalline silicates are found in 
absorption in the 10 micron complex. 

AC Her. The 10 micron complex of AC Her is one of 
the best examples of the presence of crystalline silicates. 
All crystalline silicate features are found in this source (in 
emission). This complex is well fitted by only crystalline 
silicates (paper III), but an amorphous component cannot 
be excluded. 

HD 45677: Crystalline silicate features are found in the 
10 micron complex of HD 45677. The 10 micron amor- 
phous silicate feature peaks at shorter wavelengths than 
for e.g. AC Her. 

89 Her: 89 Her is the third source showing a clear ev- 
idence of crystalline silicates in its 10 micron complex. 



Still, the weakness of the 8.3 micron feature suggests that 
the complex is dominated by the 10 micron amorphous 
silicate feature. The peak position is red shifted with re- 
spect to the ISM absorption feature (which is at 9.7 /xm), 
an indication of the presence of large grains. 

MWC 300: The 10 micron amorphous silicate absorp- 
tion feature is remarkable. There seems to be some sub- 
structure in the 10 micron complex, which might be ex- 
plained by crystalline silicate features in emission. How- 
ever, apart from the prominent 8.3 micron feature, these 
are always seen only in one scan direction and are there- 
fore not trusted. The 8.3 micron feature is very similar 
to the ones in AC Her, 89 Her and HD 45677, all sources 
which have crystalline silicate emission in the 10 micron 
complex. In MWC 300 it seems that there is amorphous 



Fig. 7. The 10 micron complex for the outflow sources, plus the mean spectra and the responsivity profile of the relevant bands 
in this wavelength region (SWS-band 2C and 3A). The gray line is the mean outflow spectrum, except in the lower right corner 
where it is the mean disk spectrum. The vertical gray lines indicate the mean peak positions of the features found, while their 
range is indicated by the errorbar on top of the plot. The thick line in the upper left corner in each panel indicates the mean 
noise level. 



silicate absorption. The E(B-V) to this star is 1.19 (Voors 
1999) which corresponds to an Ay of 3.7, assuming a nor- 
mal reddening law. This would correspond to a rg.7 — 0.25 
using the extinction law of Sandford et al. (1995). Based 
on the ISO spectra we derived rg.r = 0.13 ± 0.03. The in- 
terstellar extinction curve presented by Roche & Aitken 
(1984) would give Tg.7 — 0.2, which still above the ob- 
served Tg,r- The lower value from the ISO data suggests 
that the shape and strength of this feature is not only due 
to interstellar absorption. A natural explanation would be 
that the interstellar absorption profile is filled in by cir- 
cumstellar emission. The central wavelength of this broad 
circumstellar emission feature is likely shifted from the 



9.7 /im interstellar absorption peak, because the observed 
profile does not look like the typical interstellar absorp- 
tion profile. This shift to longer wavelengths might be ex- 
plained by large grains (see Fig 

HD 44179: The 10 micron complex of HD 44179 is 
dominated by PAH-features. No crystalline silicates seem 
to be present. A possible crystalline silicate feature might 
be seen around 10.8 /im, however at about the same po- 
sition a feature is seen in genuine PAH-sources, and at- 
tributed to a PAH cation. 

Roberts 22: The 10 micron complex is dominated by 
PAH features. The noise in the spectrum prevents us from 



Fig. 8. The Qabs for spherical amorphous ohvine grains 
(MgFeSi04; Dorschner et al. 1995) for different grain sizes: 
0.01 /im (solid line), O.l/im (dotted line), 0.3/xm (short 
dashed line), l.O/im (long dashed line), 3.0/ini (dashed 
dotted line). The curves were scaled to get an equal 
strength around 10 fim and then offset from each other. 



reaching a conclusion about the presence or absence of 
crystalline silicate features. 

4.1.2. Outflow sources 

Vy 2-2: The 10 micron complex of Vy 2-2 is dominated by 
the amorphous silicate feature, which is relatively broad. 
This feature starts at relatively blue wavelengths which 
suggests that the 8.3 micron feature is present. Because of 
the noise level (increasing with wavelength) no statements 
could be made about the presence of other crystalline sil- 
icate features. At 12.8 iim the [Ne II] line is found. 

HD 161796: As in all outflow sources, the 10 micron 
feature in HD 161796 is dominated by the 10 micron amor- 
phous silicate feature. The signal to noise level is too low 
to detect the crystalline silicate features. 

OH 26.5+0.6: The 10 micron complex is domi- 
nated by amorphous silicate absorption. Some structure 
seems present, however the absorption pattern cannot be 
matched with the features of crystalline silicates seen in 
other stars. Emission is not expected here, since at longer 
wavelengths we do see crystalline silicate features in ab- 
sorption. The origin of the substructure therefore remains 
unclear. 

HD 179821: The 10 micron complex is dominated by 
the red-shifted amorphous silicate feature. The shape and 
position resembles the feature of AFGL 4106 where it has 
been attributed to the presence of large grains (Molster et 
al. 1999b). 

AFGL 4IO6: This complex is characterized by the 10 
micron amorphous silicate feature. The sharp feature at 
11.06 /im is instrumental; no correction, as in the other 
sources, has been applied to this feature. There is no evi- 
dence for crystalline silicate features. 

NML Cyg: The 10 micron feature is dominated by 
amorphous silicate absorption. The emission wings, at 
both sides together with the location of the center of the 
absorption profile suggest that it is mainly self absorption. 



A (small) contribution from interstellar extinction cannot 
be excluded or confirmed. The substructure at 9.39, 9.54, 
10.07, 10.33, 10.76 and 11.0 /urn is due to gas-phase NH3 
(Yamamura private comm.). 

IRC+IO42O: The 10 micron feature is dominated by 
amorphous silicate emission. It peaks at a wavelength sig- 
nificantly offset from the usual 9.7 micron feature. This 
might be due to large grains, as in AFGL 4106. Fix and 
Cobb (1988) also suggested the presence of large grains 
to explain the 10 over 20 micron amorphous silicate ra- 
tio (see Fig. ^). There are indications of substructure in 
the amorphous silicate feature, which is also seen in stars 
with prominent 10 micron crystalline silicate features. In- 
terestingly enough no 11.3 micron feature (attributed to 
forsterite) is detected. This suggest that the forsterite 
abundance and/or temperature is probably very low, in 
contrast to the enstatite abundance. It should be noted 
that in the rest of the spectrum the forsterite features are 
weak or not detected. 

4.2. The 18 micron complex 

The 18 micron complexes of the disk and outflow sources 
are shown in figures ^ and |lO|. 

4.2.1. Disk sources 

IRASO9425-6O4O: The 18 micron complex of IRAS09425- 
6040 is quite similar to the mean 18 micron complex of 
the disk sources, both with regard to the positions of the 
features as well as in the strength of the features. 

NGC 6537: Because of the low dust temperature, the 
18 micron complex is less prominent than in other sources. 
Still, most prominent features are easy to recognize. The 
18.1 micron feature is relatively strong compared to the 
other features. The 18.9 and 19.5 micron features are 
nicely separated. The sharp feature at 18.88 /^m is uniden- 
tified; see the Appendix for more details on this feature. 

NGC 6302: The 18 micron complex of NGC 6302 is 
dominated by crystalline silicates. The 19.5 micron feature 
is relatively strong compared to the 18.1 micron feature. 

MWC 922: Compared to the mean disk spectrum, 
MWC 922 shows a drop around 16.5/im. It is quite similar 
to what has been found for MWC 300. In Sec |4.1.l| it has 
been argued, that the 10 micron complex of MWC 300 
is likely to be fiUedin interstellar absorption. MWC 922 
has an E(B-V)« 2 (Voors 1999), which is higher than for 
MWC 300. It is not known which fraction is due to inter- 
stellar and which is due to circumstellar reddening, but 
this high value leaves the possibility of filled in interstel- 
lar absorption open. One would expect to see evidence for 
this hypothesis in the 10 micron complex. Unfortunately 
that complex is dominated by PAH features, which will 
fill in any absorption profile and make it more difficult to 
recognize, since these features are located on the slopes of 
the silicate absorption profile, masking its presence. The 



Fig. 9. The 18 micron complex for the disk sources, together with the mean disk spectrum (gray Hne). The vertical gray lines 
indicate the mean peak positions of the features found, while their range is indicated by the errorbar on top of the plot. The 
thick line in the upper left corner in each panel indicates the mean noise level. 



absorption by crystalline silicates in the 10 micron com- 
plex suggest that circumstellar amorphous silicates will 
also be in absorption at 10 and 18 /xm. Another indirect 
argument for interstellar absorption in the 18 micron com- 
plex is the peak wavelength of the absorption feature. It 
is much bluer than the peak wavelength of the amorphous 
silicates found in emission in the disk sources. 

This is the only star where the 17 micron feature is 
clearly divided into a 16.7 and 17.0 micron feature. In 
other sources where these features are found, they are 
severely blended. A closer look at the mean disk spectrum 
indicates that this feature is a blend of 2 features, even if 
one excludes MWC 922 from this mean spectrum. It is also 
found in the two independent observations of MWC 922, 
so all these arguments make us confident about the real- 



ity of these two features. The 19.5 micron feature is the 
bluest found in our sample. There seems to be weak struc- 
ture on top of the 18.1 micron feature, which is also found 
in HD 44179. 

AC Her. The 18 micron complex is dominated by the 
15.3, 15.9 and 16.2 micron features. These features are rel- 
atively strong because of the high temperature the crys- 
talline silicates around this source attain. As in MWC 300 
and 89 Her, the 18.9 micron feature is rather blue-shifted. 
This object has a low E(B-V) due to interstellar extinc- 
tion 0.1). This low value of the interstellar E(B-V), 
the strength of the 15.9 and 16.2 micron features, and the 
fact that the spectrum returns to the same level around 15 
/im, suggests that the spectral structure around 17 /im in 
this object is caused by the 16.2 and 18.1 micron features. 



Fig. 10. The 18 micron complex for the outflow sources, plus the mean spectra, and the responsivity profile of the relevant 
bands in this wavelength region (SWS-band 3A, 3C and 3D). The gray line is the mean outflow spectrum, except in the lower 
right corner where it is the mean disk spectrum. The vertical gray lines indicate the mean peak positions of the features found, 
while their range is indicated by the errorbar on top of the plot. The thick line in the upper left corner in each panel indicates 
the mean noise level. 



Note that in this source the spectrum drops at 18 fim, 
while in MWC 300 and MWC 922 the spectrum stiU rises 
longwards of 18 ^m, indicating that in the latter sources 
the 18.1 micron feature is on a slope. 

HD 45677: The 18 micron complex of HD 45677 is very 
similar to the mean disk spectrum. It shows evidence for 
the presence of amorphous silicates. On top of the amor- 
phous silicate feature, weak crystalline silicate features are 
found. Among the disk sources, it is the object with the 
most distinct presence of amorphous silicates. The amor- 
phous silicate feature peaks at a rather red wavelength, 
which might be an indication of relatively large grains. 
However, there is only a small shift in the 10 micron com- 
plex, which suggest that the grains are relatively small. 



Therefore a compositional difference of the dust particles 
is more likely to explain the red peak position. We can 
also not exclude that there are contributions from other 
materials, such as simple oxides and/or other silicates. 

89 Her. The 18 micron complex of 89 Her shows a 

broad emission plateau, imlike most disk sources. This 
is evidence for the presence of amorphous silicates. Its 
profile is rather broad, with a prominent blue wing. On 
top of this amorphous feature, the crystalline features are 
present. The 18.0 micron feature is relatively strong. Sim- 
ilar to MWC 300 and AC Her, the 18.9 micron feature 
is severely blue-shifted (towards 18.49 /im in this source), 
which makes the 19.5 micron feature rather prominent. 
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The 16.9 micron feature is more pronounced than in most 
spectra. 

MWC 300: The 18 micron complex of MWC 300 
is characterized by (interstellar) silicate absorption to- 
gether with circumstellar silicate emission. The circum- 
stellar amorphous silicate emission profile is red-shifted 
with respect to the normal (=interstellar) amorphous sil- 
icate profile. This might indicate that large grains are 
present, a fact which can also be inferred from the high 
mm continuum flux (Henning et al. 1994). The different 
features of crystalline silicates are visible in the spectrum, 
although influenced by interstellar absorption. The clear- 
est example of a feature influenced by interstellar absorp- 
tion, is the 18.0 micron feature. In most sources it is one 
of the strongest features in the 18 micron complex, while 
in MWC 300 it is rather weak, since it has to fill up part 
of the interstellar absorption feature. The removal of the 
interstellar absorption profile is likely to reveal a promi- 
nent 18.0 micron feature in MWC 300 too. The shape 
of the 18.0 micron feature is also influenced by the blue- 
shifted 18.9 micron feature, a shift which was also found 
in AC Her and 89 Her. The 19.5 micron feature is promi- 
nently present and rather broad, as is also seen for the 
AC Her and 89 Her. In the spectrum of MWC 922 the of- 
ten blended 16.9 micron feature is split into two sharp fea- 
tures at 16.7 and 17.0 /im. In the spectrum of MWC 300 
only the sharp 17.0 micron feature is present, while the 
16.7 micron feature is missing. 

Crystalline silicates also influence the shape of this 
complex. The 18.1 and 16.2 micron features are located at 
the positions of the rise in the 'absorption' feature. This 
also leads to a change in the feature appearance and is an- 
other proof of silicate emission together with interstellar 
absorption. We exclude the possibility of no absorption at 
all and only emission features, because this would lead to 
a broad emission plateau between the 10 and 18 micron 
region, which has not been observed in other stars. Nar- 
row absorption like profiles are also found in the spectra 
of MWC 922, AC Her, IRAS09425-6040 and NGC 6302 
(see Fig.^. However, in most of these cases it is caused 
by strong 16.2 and 18.1 micron features. This is evidenced 
by the decrease of the flux (in the continuum subtracted 
spectra) on the blue side of the blend of the 15.9 and 16.2 
micron features. This drop is not seen in MWC 300 nor 
probably MWC 922. 

HD 4-in9: It cannot be excluded that part of the fea- 
tures in this complex originate from a carbon based chem- 
istry. The 18 micron complex differs in several respects 
from the mean disk spectrum. The 18.9 micron feature is 
the strongest seen in this sample. The spectrum around 
15 /um is influenced by CO2 absorption at 13.9, 14.9, 15.3 
and 16.2 /um. Without these absorptions the feature at 
15 /zm would be very similar to the one found in 89 Her. 
The 15.9 and 16.2 micron features are rather prominent, 
although not as strong as in AC Her. The 16.5 micron 
feature, probably a PAH-feature (Van Kerckhoven et al. 



2000), is present in both the AOTOl and the AOT06 spec- 
trum and seems therefore rather secure. The only other 
star where this feature has been detected is NGC 6537. 
Also, longwards of 16.5 /im these two stars look quite sim- 
ilar, with similar features at 16.9, 18.0 and 18.9 /xm. The 3 
features found around 18 /um (at 17.8, 18.0 and 18.2 ^m) 
are seen in both the AOTOl and the AOT06 spectra as 
well as in the up and the down scans and there is not 
much doubt about their reality. The 19.5 micron feature 
is rather weak in HD 44179. 

Roberts 22: The 18 micron complex is rather noisy, pre- 
venting the measurement of the individual components. 
The gentle slope at the blue side of the amorphous silicate 
feature suggests the presence of some crystalline silicates 
around 16 fMa. Other crystalline silicate features might be 
present but the noise prevents their detection. 

4.2.2. Outflow sources 

Vy 2-2: Unfortunately no interesting features could be 
confirmed in this wavelength range, because of the low 
signal to noise ratio. The emission line at 15.56 /im is 
[Ne HI]. 

HD 161796: The 18 micron complex of HD 161796 
is dominated by the 18.1 micron feature. However, one 
should be careful with the interpretation of this feature 
(see Appendix for more details). The 19.5 micron feature 
is relatively strong. 

OH 26.5+0.6: The 18 micron complex of OH 26.5-^0.6 
is an absorption spectrum. Both the amorphous and the 
crystalline silicates are seen in absorption. We can rec- 
ognize most of the features which are normally seen in 
emission. All absorption features are slightly red-shifted, 
we speculate that this red-shift is due to a temperature 
effect (Bowey et al. 2000). It is clear from the absorption 
profile that crystalline silicates only play a minor role in 
the spectrum (at these wavelengths). 

HD 179821: This complex lies on a very steep slope. It 
shows the characteristic outflow source spectrum. It much 
resembles AFGL 4106 and IRC-l-10420, two other massive 
stars. It is dominated by the 18 micron amorphous silicate 
feature. On top of this feature, weak crystalline silicates 
can be detected. 

AFGL 4IO6: The 18 micron complex of AFGL 4106 is 
very similar to the mean outflow spectrum, a broad amor- 
phous silicate feature with weak crystalline silicate fea- 
tures on top. No differences between the mean spectrum 
and AFGL 4106 were found. 

NML Cyg: The 18 micron amorphous silicate absorp- 
tion is less prominent than for similar high mass star 
IRC-l-10420. This might indicate that the 18 micron amor- 
phous silicate band is just in the regime between self- 
absorption and emission. Whatever the reason, it makes 
the crystalline silicates more apparent. Most crystalline 
silicate features are clearly seen. Still, the 19.5 micron fea- 
ture seems absent. 
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IRC+IO42O: The 18 micron spectrum is very similar 
to the mean outflow spectrum. A strong amorphous com- 
ponent is seen with some weak crystalline silicate features 
on top of it. 

4-.3. The 23 micron complex 

The 23 micron complex spectra of the disk and outflow 
sources are presented in figures |ll| and 

4.3.1. Disk sources 

IRASO9425-6O4O: As for the 18 micron complex, the 23 
micron complex in this source is also quite similar to the 
mean 23 micron complex for the disk sources. There are 
some differences: The slope to the short wavelength side 
continues to shorter wavelengths than in the mean disk 
spectrum and the plateau does not show the step-like 
structure seen in the mean (disk) spectrum. The 24.5 mi- 
cron feature is relatively strong. 

NGC 6537: The 23 micron complex of NGC 6537 is 
characterized by a very strong 23.7 micron feature and 
relatively weak 23.0 micron feature. Although the 24.5 mi- 
cron feature is severely influenced by the [Ne V] emission 
line, the step-hke structure is visible in the plateau. 

NGC 6302: The main differences between the mean 23 
micron complex and the 23 micron complex of NGC 6302, 
is found in the strength of the 23.0 and 25.0 micron fea- 
tures. The 23.0 micron feature is slightly weaker than aver- 
age while the 25.0 micron feature is significantly stronger. 

MWC 922: The most striking aspects of the 23 micron 
complex of MWC 922 is the strong 23.9 micron feature 
and thus the relative weakness of the 23.7 micron feature. 
Also the 23.1 micron feature is relatively strong. However, 
there is no correlation between these two features; e.g., 
compare MWC 922 with AC Her. 

AC Her. The 23 micron complex differs in several re- 
spects from the mean spectrum. The feature at 22.69 /im 
has a remarkable shape: it is sharp and rather prominent. 
It appears very similar to a feature in HD 44179, but there 
it is weaker. However, there are some questions about the 
reality in AC Her (see the Appendix) . If this feature can 
be attributed to the 23.0 micron feature it is the most 
blue-shifted, which implies that also the 23.7 micron fea- 
ture is significantly blue-shifted. Alike to MWC 922, there 
is a 23.89 micron feature present. The 24.5 micron feature 
in the spectrum of AC Her is not as prominent as in the 
mean spectrum. 

HD 45677: The 23 micron complex of HD 45677 is 
rather similar to the mean 23 micron disk complex. Still 
there are two subtle differences. The 25.0 micron feature 
is weaker in HD 45677 than in the mean spectrum, while 
at the blue side of this complex the 23.0 micron feature 
is stronger than in the mean spectra. This is exactly the 
opposite of what is found in NGC 6302. 



89 Her. The 23 micron complex of 89 Her is rather 
smooth and few individual structures can be recognized. 
The broadness of this complex is similar to the complexes, 
including the plateau, seen in other sources. The width of 
the features seems larger than in the average disk spec- 
trum and this may cause severe blends. Still, better signal- 
to-noise will probably unveil similar structures to those 
found in the other sources, 

MWC 300: The 23 micron complex attracts attention 
because of the absence of the red shifted 23.0 micron fea- 
ture. Two of the components of the 25 micron plateau are 
quite clear in this object. 

HD 441^9: Where the 18 micron complex looks quite 
similar to NGC 6537, the 23 micron complex is rather 
different. It has a rather prominent 23.0 micron feature 
and the slope at the red side of the 23.7 micron feature is 
very gentle. There are indications for the step-function of 
the plateau but the noise level prevents strong conclusions. 
The 23.0 micron feature is rather prominent and blue- 
shifted with respect to the mean spectrum. 

Roberts 22: The 23 micron complex is dominated by 
noise and by fringes. The strong structures around 23 /^m 
and 25 /im are due to constructive interference of the 
fringes in the up and down scan, and the relatively fiat 
structure around 24 /tm is due to destructive interference 
of the fringes. The width of the whole structure is quite 
similar to the mean spectrum width. This suggests that 
the complex is present but not much can be said about 
the structure of this complex. 

4.3.2. Outfiow sources 

Vy 2-2: The spectrum in this wavelength range also suffers 
from a lot of noise. The comparison with the mean outflow 
spectrum suggests the presence of the 23.0 and 23.7 micron 
blend. 

HD 161796: Apart from the 23.1 and the 23.7 micron 
features, the most striking feature is the 25.0 micron fea- 
ture. It is rather blue-shifted with respect to the 25.0 mi- 
cron features in other sources. Its presence is also con- 
firmed in the rev071 dataset, although the feature there 
seems less blue-shifted and more in line with the other 
25.0 micron features. 

OH 26.5+0.6: As in the case of the 18 micron complex, 
the crystalline features appear as weak absorption features 
on top of the wing of the amorphous silicate absorption. 
The blend of the 23.0 and 23.7 micron features is rather 
narrow compared to the mean spectrum. The 23.89 micron 
feature seems present in absorption. 

HD 179821: As for all outflow sources, the 25 micron 
plateau is also very weak in this source. The feature at 23.0 
micron feature is surprisingly strong and the 23.7 micron 
feature rather broad. The sharp rise of the 23.0 micron 
feature, which is typical for the outflow sources, is evident 
here. 



Fig. 11. The 23 micron complex for the disk sources, together with the mean disk spectrum (gray hue). The vertical gray lines 
indicate the mean peak positions of the features found, while their range is indicated by the errorbar on top of the plot. The 
thick line in the upper left corner in each panel indicates the mean noise level. 



AFGL 4IO6: The 23 micron complex is somewhat re- 
markable. The 25 micron plateau is, as for most outflow 
sources, rather weak and ends already before 25 /im. The 
features around 23 /im produce a very flat plateau. Some 
structure can be detected, but together with IRAS09425- 
6040 it shows the flattest structure in our whole sample. 

NML Cyg: The 23 micron complex is interesting in two 
ways. The 25 micron plateau is almost absent. However, 
the blend of the other features is the broadest seen in our 
sample. In this respect the spectrum is somewhat similar 
to the 23 micron complex of IRC+10420. 

IRC+IO42O: The 23 micron complex is weak with re- 
spect to the continuum and therefore affected by noise. 
Still we could identify the different features. The complex 



is also quite broad and alike to what is seen for AFGL 4106 
and NML Cyg. 

4-. 4- The 28 micron complex 

The 28 micron complex spectra of the disk and the outflow 
sources are plotted in figures ^ and 

4.4.1. Disk Sources 

IRASO9425-6O4O: The strength of its 27.6 micron feature 
is an intriguing aspect of the 28 micron complex. It is also 
very broad and blends with the 26.9 and 28.2 micron fea- 
tures. The 30.6 micron feature is relatively strong. There 



Fig. 12. The 23 micron complex for the outflow sources, plus the mean spectra, and the responsivity proflle of SWS-band 3D. 
The gray line is the mean outflow spectrum, except in the lower right corner where it is the mean disk spectrum. The vertical 
gray lines indicate the mean peak positions of the features found, while their range is indicated by the errorbar on top of the 
plot. The thick line in the upper left corner in each panel indicates the mean noise level. 



is a hint of the 31.2 micron feature in the rev084 data. 
However, this cannot be confirmed by the rev254 data 
because of the noise level in this area. The 29.6 micron 
feature peaks at the very blue side of the feature in the 
mean disk spectrum. 

NGC 6537: One of the most interesting aspects of the 
28 micron complex in NGC 6537 is the absence of the 
31.2 micron feature. The 29.6 and 30.6 micron features 
are both very strong relative to the 27.6 and 28.2 micron 
features. 

NGC 6302: The 29.6 and 30.6 micron features are very 
strong compared to other sources, but very similar to 
NGC 6537. In contrast to NGC 6537, the 31 micron fea- 
ture is prominent in this source. The strength ratio of the 



27.6 micron feature to the 28.2 micron feature is rather 
low compared to the average. 

MWC 922: The 27.6 and 28.1 micron features are rel- 
atively strong in this star. Compared to other spectra, the 
29.6 micron feature 'misses' intensity at the long wave- 
length side. If the 29.6 micron feature is indeed a blend 
of two features, then the reddest feature is severely de- 
pressed in MWC 922. The 30.6 and 31.2 micron features 
are clearly present and relatively strong. 

AC Her. The 28 micron complex is rather noisy, espe- 
cially longwards of 29 micron, where band 4 starts, but 
all features seem present. The 27.6 and 28.2 micron fea- 
ture arc still well visible and rather strong. The 29.6 and 
30.6 micron features are hardly detectable above the noise. 



Fig. 13. The 28 micron complex for the disk sources, together with the mean disk spectrum (gray hne). The vertical gray lines 
indicate the mean peak positions of the features found, while their range is indicated by the errorbar on top of the plot. The 
thick line in the upper left corner in each panel indicates the mean noise level. 



The sharp and strong peak at 31.0 is questionable. 
Still, there are indications of an underlying 31.2 micron 
feature. 

HD 45677: The 28 micron complex of HD 45677 is 
interesting due to the strong 27.6 and 28.2 micron features. 
The ratio between the flux of these two features and the 
29.6 micron feature is the highest found in our sample. 
The 30.6 micron feature is also weak, therefore the dust 
species that causes the prominent 29.6 and 30.6 micron 
features in NGC 6302 and NGC 6537 must be only a minor 
component in HD 45677. 

89 Her. The 28 micron complex is dominated by the 
rise of the broad 33 micron band. On top of this feature are 
the 27.6 and 28.2 micron features. We have possibly foimd 
a broad structure around 29.7 /xm and no indications for 



a 30.6 micron feature. The 33 micron band together with 
the low S/N for Band 4 make this a problematic part of 
the spectrum. 

MWC 300: MWC 300 is again a source where the 27.6 
and 28.2 micron features are stronger than the 29.6 and 
30.6 micron features. The 28.8 micron feature is present in 
this source, but the sharp peak (at 28.7 /um) on top of the 
28.8 micron feature has significant unccrtaintic^s (sec Ap- 
pendix). Neglecting this peak, the remaining feature lines 
up very well with the weak feature found in NGC 6537. 

HD 441'^9: The 28 micron complex is, like that of 
89 Her, dominated by the onset of the broad 33 micron 
band. On top of the blue rise of this broad band, the 27.6 
and 28.2 micron features are rather prominent (in con- 
trast, for example, to NGC 6302). The 29.6 micron feature 



Fig. 14. The 28 micron complex for the outflow sources,, plus the mean spectra, and the responsivity profile of the relevant 
bands in this wavelength region (SWS-bamd 3D, 3E and 4). The gray line is the mean outflow spectrum, except in the lower 
right corner where it is the mean disk spectrum. The vertical gray lines indicate the mean peak positions of the features found, 
while their range is indicated by the errorbar on top of the plot. The thick line in the upper left corner in each panel indicates 
the mean noise level. 



seems less pronounced and also the 30.6 micron feature is 
just above the noise level. The 31.2 micron feature is not 
detected. 

Roberts 22: Because of the flux jumps around band 3E, 
it was difficult to subtract a reliable continuum for band 
3E. Wc have tried to maintain the relative slope after a 
jump in the flux. In this way we were able to subtract the 
continuum, but the relative strength of the features in the 
different bands (3D, 3E and 4) in this wavelength range is 
severly influenced by these flux jumps. We measured the 
bands in this wavelength range using a local continuum. 
Wc foimd the 27.6, 28.2 and 28.8 micron features, all three 
of which were seen in both the rev254 and the rev084 



spectra. The 28.8 micron band is interesting since it is the 
most red-shifted one from the four stars where this band 
is found. 

4.4.2. Outflow Sources 

Vy 2-2: The beginning of band 4 is less effected by noise. 
Band 4 has a larger aperture and is therefore less influ- 
enced by the mispointing. The 29.6 and 30.6 micron fea- 
tures can be identified. 

HD 161796: The noise is rather severe in this part of 

the spectrum. Still, some features can be identified. The 
blend of 27.6 and 28.2 micron features is present. The 
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29.6 micron feature is detected but has an extension at 
the short wavelength side, in both the up and down scan, 
which is not seen in the rev071 data. This leaves some 
uncertainty on the reality of this extension. 

OH 26.5+0.6: This is the complex where the absorp- 
tion spectrum of this source goes into an emission spec- 
trum. The 29.6 micron feature, the first feature in this star 
which appears in emission, is rather blue-shifted. No 30.6 
or 31.2 micron feature is found. 

HD 179821: The 27.6 and 28.2 micron features are 
relatively weak with respect to the 29.6 micron feature. 
Whether this is a typical outflow source characteristic 
is difficult to say, because the mean outflow spectrum is 
rather influenced by this spectrum (see paper II). 

AFGL 4IO6: because of badly corresponding up and 
down scans we have removed the first part of band 4 (29 
to 30.5 /xm). Still, the 27.6 and 28.2 micron features stand 
out very prominently, with sharp drops on either side of 
the blend. 

NML Cyg: The end of band 3D and the beginning of 
Band 3E suffer from spectral leakage from shorter wave- 
lengths, therefore no features are measured in this wave- 
length range. The 30.6 micron feature is prominent in this 
spectrum, while, on the other hand, there is no indication 
of the 31.2 micron feature. 

IRC+IO42O: Part of the spectrum has been removed 
due to deviating up and down scans. The 27.6 and 28.2 mi- 
cron features are clearly present, but rather blue-shifted. 
As in many outflow sources, there is no evidence for a 30.6 
or 31.2 micron feature. 

4.5. The 33 micron complex 

The 33 micron complexes of the disk and the outflow 
sources are plotted in figures ^ and |l^. 

4.5.1. Disk Sources 

IRASO9425-6O4O: The 33 micron complex follows the 
mean spectrum very well. There are two major discrep- 
ancies, the 32.2 micron feature is stronger than in the 
average spectrum and the plateau drops more steeply at 
the end, which might be an indication of the 36.5 micron 
forsterite feature. It should be noted that the peak to con- 
tinuum ratio of the 33.6 feature is the highest found so far 
in the ISO database. 

NGC 6537: The 33 micron complex of NGC 6537 is 
characterized by sharp features. The 33.6 micron feature 
is very prominent The peak caused by the 32.8 micron 
feature is rather sharp, and much weaker than the 33.6 
micron peak. The plateau ends at almost 37 /xm due to the 
presence of a 36.5 micron feature. The 34.1 micron feature 
is seen as a small plateau on the 33.6 micron feature. 

NGC 6302: In general, the complex in this source is 
quite similar to the mean complex structure. However, the 
32.8 micron feature is weaker than average in this source. 



MWC 922: The 33 micron complex of MWC 922 is 
characterized by a prominent plateau. The weak structure 
(at about 33.5 /xm) on top of the 33.6 micron structure is 
likely the [S III] emission line. The 32.8 micron feature is 
slightly stronger than average. The feature at 34.8 /im is 
the [Si II] line. 

AC Her: It should be noted that the band 4 data is 
of low quality, although a lot of spurious points (glitches 
and glitch-tails) have already been removed by hand. As 
in the 23 micron complex there is some doubt about one 
of the most prominent features, the one at 32.6 /im (see 
for more details the Appendix). We note that the 33 and 
23 micron complexes have a similar appearance. In the 23 
micron spectrum of this source, a sharp feature has also 
been found at the short wavelength side. 

HD 45677: The 33 micron complex of HD 45677 is 
that of a typical disk source. It has a rather strong 33.6 
micron feature and a relatively weak 32.8 micron feature. 
The main difference is found in the strength of the 35.9 
micron feature, which is much stronger in this source than 
in the mean spectrum. As for NGC 6537, the red side of 
the 33.6 micron feature shows evidence for the 34.1 micron 
feature. The plateau extends beyond 37 /im. 

89 Her. Interpretation of the 33 micron complex suffers 
from the low ffux levels. The three main features (the 32.8 
and 33.6 micron features and the plateau) are present. Any 
substructure is difficult to quantify due to the low S/N. 
The 32.8 and 33.6 micron features are rather blue-shifted. 

MWC 300: The overall shape is very similar to the 
mean spectrum. The 33.6 and 34.1 micron bands are 
clearly split in this source, and the structure nicely corre- 
lates with the structure seen in AFGL 4106. The plateau 
extends to almost 37 /im and the sharp feature at 34.81 
/im is the [Si II] line. The 32.8 micron feature is weak but 
clearly present. 

HD 441^^9: The 33 micron complex can be described 
as a very broad feature peaking at 33.6 /im, a sharp 32.8 
micron feature and a plateau extending to wavelengths 
larger than 37 /tm. This very broad feature is a blend of 
the 33.6 and 34.0 micron features. In the 35 micron plateau 
the different components can be identified. 

Roberts 22: The general shape is roughly equal to the 
mean disk spectrum. We masked the unreliable narrow 
(0.1 iim) peak at 33.3 /tm (see Appendix for more de- 
tails), but still the 33.6 micron feature peaks at rather 
blue wavelengths like 89 Her, which shows a similar com- 
plex. The 36.0 micron feature, which is mostly found in 
outflow sources, is rather strong. It was already stated be- 
fore that this source shows a lot of outflow characteristics. 
On the other hand, the strength ratio of the 32.8 micron 
feature to the 33.6 micron feature is more like that of the 
disk sources. 



Fig. 15. The 33 micron complex for the disk sources, together with the mean disk spectrum (gray hne). The vertical gray lines 
indicate the mean peak positions of the features found, while their range is indicated by the errorbar on top of the plot. The 
thick line in the upper left corner in each panel indicates the mean noise level. 



4.5.2. Outflow Sources 

Vy 2-2: The 33 micron complex suffers also from a low 
signal-to-noisc ratio, but the general structure can be 
distinguished. The ratio of the 32.8 micron feature to 
the 33.6 micron feature is rather low in comparison with 
the rest of the outflow sources. The plateau is relatively 
strong. 

HD 161796: The 33 micron complex of HD 161796 is 
more like the disk sources than the outflow sources. The 
33.6 micron feature is stronger than the 32.8 micron fea- 
ture. Also the gentle slope at the red side of this last blend 
resembles the disk sources better than the outflow sources. 
On the other hand, the 36.0 micron feature is more typical 
of the outflow sources. 



OH 26.5+0.6: The 33 micron complex of OH 26.5+0.6 
is rather similar to that of HD 161796; especially the 
broadness of the whole structure. The main difference be- 
tween these two stars is found in the 35 micron plateau. 
The 36.0 micron feature seems present in both sources, 
but there is some extra absorption around 35.0 /xm in 
OH 26.5+0.6. 

HD 179821: The 33 micron complex is characterized by 
the two relatively sharp features at 32.8 and 33.6 micron 
and, as in the 23 micron complex, a weak plateau. The 33.6 
micron feature is one of the sharpest in our sample. We 
note that there seems to be no evidence for the 36.5 micron 
feature in both datasets. The strength ratio of the 32.8 to 
33.6 micron features suggests that this star is massive and 
likely a post-red Supergiant. 



Fig. 16. The 33 micron complex for the outflow sources,, plus the mean spectra and the responsivity proflle of SWS-band 4. 
The gray line is the mean outflow spectrum, except in the lower right corner where it is the mean disk spectrum. The vertical 
gray lines indicate the mean peak positions of the features found, while their range is indicated by the errorbar on top of the 
plot. The thick line in the upper left corner in each panel indicates the mean noise level. 



AFGL 4 106: At first sight, tlie 33.6 micron feature 
looks much stronger than the 32.8 micron feature, in real- 
ity they do not differ so much, because of the presence of 
a broad 35 micron band in the spectrum of AFGL 4106. 
This is in accordance with most other outflow sources, 
where equal strengths for these two features arc found. 
The broad underlying 35 micron band might be similar to 
the one seen in 89 Her and HD 44179. While it seems to 
start at a longer wavelength than in the other two sources, 
this might be an artifact because the data around 30 /xm 
have been ignored in this analysis. 

NML Cyg: The spurious 32.97 micron feature gives the 

33 micron complex its unique appearance. Apart from this 
detector respons structure the 32.97 micron feature looks 



normal. The 34.1 micron feature is very clear. Because 
the 33.6 and 34.1 micron features are found in both the 
up and down scans of the AOT06 and the AOTOl data, we 
are confident about their reality. This makes this source 
one of the stars where they are most clearly separated. 
At 34.6 /im, the OH mascr pump line can be found in 
absorption as was already reported by Justtanont et al. 
(1996b). 

IRC+IO42O: The 33 micron complex is dominated by 
the strong 32.97 micron feature, which is predominantly 
seen in high flux sources, both in the carbon-rich as well 
as in the oxygen-rich dust environments. Therefore we as- 
sume that it is an artificial feature caused by the detector 
responsivity for high fiuxes. Underneath this feature, the 
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usual 33 micron complex is visible, although at a rather 
low level. The 34.1 micron feature seems rather strong 
with respect to the 33.6 micron feature, although this last 
feature is also influenced by the presence of the 32.97 mi- 
cron feature. Alike to the features found in the 28 micron 
complex, these features are also blue-shifted with respect 
to the mean spectrum. In the 35 micron plateau, in which 
the different components can be detected, the 34.6 fj,m OH 
maser pumping line is nicely detected in absorption (see 
Sylvester et al. 1997). 

4-. 6. The 40 micron complex 

The 40 micron complexes of the disk and outflow sources 
are plotted in figures |l| and |l|. 

4.6.1. Disk Sources 

IRASO9425-6O4O: The 40 micron complex is similar to the 
mean. The most prominent difference is in the strength of 
the 43.0 micron feature, which is also very blue shifted. 
This suggests a significant (cold) crystalline water-ice con- 
tribution. This would be strange for a disk source and it 
is very unfortunate that this claim cannot be confirmed 
by a spectrum of the 60 micron complex. 

Another difference is the long wavelength side of the 
40.5 micron feature. The short wavelength side of this fea- 
ture fits the mean spectrum very well, and peaks at a 
rather blue wavelength. This feature might be a blend 
(see Paper II). While noisy, this might imply that in 
IRAS09425-6040 the short wavelength component is dom- 
inant over the long wavelength component. 

NGC 6537: The 40 micron complex is rather similar 
to the mean disk complex. The main difference is that the 
43.8 micron feature is somewhat weaker than on average. 
The 41.7 micron feature is rather pronounced in this star. 

NGC 6302: The 40 micron complex is also rather simi- 
lar to the mean disk spectrum. The features in NGC 6302 
seem a little sharper than in the mean spectrum. Because 
of the sharpness of the features, the enstatite and crys- 
talline water ice peaks around 43 /im are nicely separated. 

MWC 922: The 40 micron spectrum of MWC 922 is 
rather similar to the mean disk spectrum. All features 
in the mean spectrum are also seen in MWC 922, and 
the feature strength ratios are comparable. Only the 44.7 
micron feature seems somewhat stronger, but this might 
be due to the continuum subtraction. 

A C Her: The noise in this part of the spectrum makes 
it very difficult to say more about this feature than that 
the 40.5 micron feature and the 43.8 micron feature are 
present. There seems some evidence for the 43.0 micron 
crystalline water ice feature, which is supported by the 
LWS data, showing evidence for the 60 micron feature. 
The 44.7 micron feature is, also visible in the LWS spec- 
trum. 



HD 45677: The 40 micron complex does show the usual 
features, i.e. the 40.5, 41.8, 43.0, 43.8 and 44.7 micron fea- 
tures. The 43.0 and 43.8 micron feature are rather narrow. 
The 43.0 micron feature is relatively weak. Unfortunately, 
there is no LWS data to check for the presence and shape 
of the 60 micron feature. 

89 Her. The 40 micron complex is completely domi- 
nated by noise and therefore no conclusions can be drawn. 

MWC 300: The 40 micron complex of MWC 300 suffers 
from the low fiux level in this part of the spectrum. Still 
we can disentangle the 40.5 micron feature, the blend of 
the 43.0 and 43.8 micron features and the 41.8 micron 
feature. 

HD 44179: The 40 micron complex of HD 44179 is 
interesting because it shows a prominent 40.5 micron fea- 
ture, but the 43.0 and 43.8 micron features are much less 
pronounced. The steep slope at the red side suggests that 
they are present, but part of the spectrum between the 
40.5 and 43.0 micron features must have been filled in. 
This might be due to an exceptionally strong 41.8 micron 
feature, but the width and wavelength of this feature as 
measured in the other stars indicate that it is not suffi- 
cient to solely account for the closing of the gap between 
the 40.5 and 43.0 micron features. 

Roberts 22: The 40 micron complex is dominated by 
the 43.0 micron feature, indicating that there is a lot of 
crystalline water ice. In this respect the spectrum looks 
somewhat similar to that of HD 161796, although the rel- 
ative strength of the 40.5 micron feature in Roberts 22 is 
much stronger than in HD 161796. 

4.6.2. Outfiow Sources 

Vy 2-2: The low signal to noise ratio in this part of 
the spectrum prevents us from reaching firm conclusions 
about the structure in this complex. It seems that the 40.5 
micron feature and the 43.0 + 43.8 micron blend can be 
distinguished. 

HD 161796: The 40 micron complex is completely 
dominated by a very strong 43.0 micron (crystalline wa- 
ter ice) feature. There might be a red-shifted 43.8 micron 
feature in the strong 43.0 feature, but this is difficult to 
quantify. Most of the features found in other 40 micron 
complexes are present in this complex. There is an indi- 
cation that the 40.5 micron feature is in fact a blend of 
2 separate bands. Closer inspection of the spectrum sug- 
gests that there is a broad band underlying the strong 
crystalline H2O band. A suitable candidate for this fea- 
ture might be amorphous H2O ice. This helps the fitting 
of the 40.5 micron band (Paper HI). 

OH 26.5+0.6: The spectrum of OH 26.5+0.6 looks 
very similar to the mean outfiow spectrum. Where the 
spectra of OH 26.5-1-0.6 and HD 161796 were very simi- 
lar for the 33 micron complex, they look quite different 
in the 40 micron complex. The strength ratio of the 40.5 
and 43.0 micron features is dramatically different. Still, as 



Fig. 17. The 40 micron complex for the disk sources, together with the mean disk spectrum (gray hne). The vertical gray lines 
indicate the mean peak positions of the features found, while their range is indicated by the errorbar on top of the plot. The 
thick line in the upper left corner in each panel indicates the mean noise level. 



for all outflow sources, the 43.0 micron feature is stronger 
than the 43.8 micron feature. The 44.7 micron feature is 
also weakly present in this spectrum. 

HD 179821: Like HD 161796, the spectrum of 
HD 179821 is dominated by the 43.0 micron crystalline wa- 
ter ice feature. However, the 40.5 micron feature is much 
stronger, rather red-shifted, and shows a very gentle slope 
on the blue side. The 43.8 enstatite feature produces some 
recognizable structure. 

AFGL 4IO6: The 40.5 micron feature looks strong, 
but again this is partly due to the broad feature peak- 
ing around 35 /um. As for HD 179821, the 40.5 micron 
feature has a gentle slope; a characteristic feature of out- 
flow sources. The 43.8 micron feature is weak and almost 



completely overwhelmed by the 43.0 micron water ice fea- 
ture. 

NML Cyg: The general shape of the complex is quite 
similar to what has been found for other sources. The 
40.5 micron feature is relatively strong with respect to 
the 43.0 + 43.8 micron blend. However this is seen in 
more sources, with HD 44179 as the most extreme case. 
The gas-phase rotational water lines, seen in emission, are 
most striking in this wavelength region (Justtanont et al. 
1996b). The modeling of these lines is outside the scope of 
this pap(T. Besides gas phase water, crystalline water ice 
is also present, as evidenced by the presence of the 43.0 
micron feature. 

IRC+IO42O: The 40 micron complex suffers more from 
the noise than the other complexes. Still the main features 



Fig. 18. The 40 micron complex for the outflow sources, plus the mean spectra and the responsivity profile of SWS-band 4. 
The gray line is the mean outflow spectrum, except in the lower right corner where it is the mean disk spectrum. The vertical 
gray lines indicate the mean peak positions of the features found, while their range is indicated by the errorbar on top of the 
plot. The thick line in the upper left corner in each panel indicates the mean noise level. 



can be distinguished. The 42.8 micron feature is difficuh 
to disentangle, but there seems hardly any 43.6 micron 
feature. 

^.7. The 60 micron complex 

The 60 micron complexes of the disk sources and the out- 
flow sources are plotted in figures ^ and |2^. 

4. 7.4. Disk Sources 

NGC 6537: The 60 micron complex is rather smooth and 
peaks at about 65 /xm. The 69 micron forsterite feature 
is clearly present. The gentle rise at the blue side of this 



feature indicates that crystalline water ice is present, as 
was already indicated by the 43.0 micron feature in the 40 
micron complex. Another dust component is necessary to 
explain the peak position and the red side of this feature, 
and diopside might be a good candidate for this (Koike 
et al. 2000). The forbidden emission lines were removed 
by hand, but in the original spectrum the following lines 
were found: [O III] at 51.79 fim, [N III] at 57.31 ^lm and 
[O I] at 63.16 ^m. It is interesting to note that all lines 
are blue shifted (in both scan directions). We probably see 
here the fast outflow from the star. Based on the shifts, the 
outflow velocity in our direction is 135±30 km/sec. This is 
comparable to the velocity found by Cuesta et al. (1995). 
They came up with a model, based on [N II] profiles at 



Fig. 19. The 60 micron complex for the disk sources, together with the mean disk spectrum (gray hne). The vertical gray lines 
indicate the mean peak positions of the features found, while their range is indicated by the errorbar on top of the plot. The 
thick line in the upper left corner in each panel indicates the mean noise level. 



6584 A, which predict a velocity of 180 ± 40 km/sec in the 
direction of the observer (400 ± 100 km/sec total outflow 
velocity). 

NGC 6302: The 60 micron complex is very similar to 
the mean disk spectrum and lines up pretty well with those 
of NGC 6537 and MWC 922. Compared to MWC 922, 
NGC 6302 shows a little more intensity on the blue side, 
which probably indicates that there is a larger contribu- 
tion from crystalline water ice as is expected from the 40 
micron complex. The 69.0 micron forsterite feature is very 
prominent and quite narrow. 

MWC 922: The 60 micron complex of MWC 922 peaks 
at very long wavelengths, implying that crystalline water 
ice does not play an important role at these wavelengths. 
This is in contrast with the 40 micron complex, where a 



band at 43.0 micron has been found. In Paper III we will 
show that the 40 micron complex of this source can nicely 
be fitted with only enstatite, confirming the absence of 
crystalline water ice. 

The 69 micron feature is very prominent and seems 
broadened at its red side, which might be an indication 
for another dust component. The rather sharp 65 micron 
feature seems flanked by another feature at 66.5 /um. 

HD 441'^9: The 60 micron complex is rather flat and 

does not show much structure, apart from the prominent 
69 micron feature and the [O I] line at 63.18 jiui. Because 
of this broad flat complex it is likely that there is only 
a small amount of crystalline water ice together with an 
extra dust component radiating at these wavelengths. 



Fig. 20. The 60 micron complex for the outflow sources, plus the mean spectra and the responsivity profiles of LWS band 1 
(solid line), band 2 (dashed line), band 3 (solid line) and band 4 (dashed line). The gray line is the mean outflow spectrum, 
except in the lower right corner where it is the mean disk spectrum. The vertical gray lines indicate the mean peak positions of 
the features found, while their range is indicated by the errorbar on top of the plot. The thick line in the upper left corner in 
each panel indicates the mean noise level. 



Roberts 22: While most disk sources do peak at a wave- 
length around 65 /im, Roberts 22 seems to peak around 
60 /zm. This is very much alike to the outflow sources and 
HD 161796 in particular. The 63.18 fj,m [O I] line is also 
visible. In correspondence with the outflow sources, there 
is hardly any 69 micron feature. 



4.7.2. Outflow Sources 

Vy 2-2: The 60 micron complex is typical of an outflow 
spectrum. It peaks around 60-62 /im and the 69 micron 
feature is absent. 



HD 161796: The 60 micron spectrum is characterized 

by the crystalline water ice structure. It peaks at w 60/im, 
which is much bluer than the disk sources, but in better 
agreement with the expectation that there is crystalline 
water ice aroiind this star. There is a hint of a 69 mi- 
cron feature. We have measured it, but better data seem 
necessary for this feature. 

OH 26.5+0.6: The 60 micron complex has 2 peaks, one 
at about 60 /xm, probably due to crystalline water ice, and 
one at 65 jxm, more like MWC 922. The 69 micron feature 
is clearly present, which is somewhat unexpected for an 
outflow source. 
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HD 179821: The 60 micron complex peaks at rather 
blue wavelengths, indicating the presence of a lot of crys- 
talline water ice; as was already concluded from the 40 mi- 
cron complex. 

AFGL 4 106: The 60 micron complex peaks around 
60 micron. This is typical for sources which are dominated 
by crystalline water ice, and in agreement with the 40 mi- 
cron complex spectrum. 

NML Cyg: The 60 micron feature is very similar to 
the mean outflow spectrum, indicating that most of this 
feature is formed by emission from crystalline water ice. 

IRC+IO42O: The 60 micron complex peaks at very 
short wavelength, indicating that there is crystalline water 
ice and a relatively low abundance of crystalline silicates, 
as was also found in the other complexes. 

4-8. Remaining Bands 

We chose not to show the individual remaining bands of 
our programme stars. However, in Paper II we show the 
average band profiles of the disk and outflow sources of 
the most prominent remaining bands, near 14, 21, 26, 48 
and 90 /zm. 

4.8.1. Disk Sources 

IRASO9425-6O4O: The (instrumental) features around 
14 fj.m are weak or absent in IRAS09425-6040. The 20.7 
and 21.5 micron features are clearly distinguishable in the 
spectrum. The 26.1 micron feature is also rather promi- 
nent. We unfortunately have no data for features at wave- 
lengths longer than 45 jjui. 

NGC 6537: As in IRAS09425-6040, the (instrumental) 
features around 14 /im arc weak or absent in NGC 6537. 
The 20.7 and 21.5 micron feature are present. The 26.1 
micron feature is one of the reddest in our sample. The 
47.7 and 48.6 micron features are separately identified. 
Finally there is clear evidence for the presence of a 91 
micron feature, with the [O III] line at 88.35 fim. 

NGC 6302: NGC 6302 has very nice 20.7 and 21.5 mi- 
cron features. The 26.1 micron feature is obscured by the 
presence of a strong emission line ([O IV]) and its pres- 
ence cannot be confirmed at the moment. The 48 micron 
features is split up into its separate components, i.e. the 
47.7 and 48.6 micron features. The 91 micron feature is 
present and even rather prominent. 

MWC 922: The 3 (instrumental) features around 
14 micron are all detected in this source. The 20.7 and 
21.5 micron arc clearly present and arc only slightly less 
prominent than in IRAS09425-6040. The 26.1 micron fea- 
ture is also rather prominently present. 

AC Her: The 20.5 and 21.5 features are more blended 
than in other stars and the final result looks somewhat 
similar to the responsivity profile. This might point to a 
problem in the dark-current subtraction. However, there 
are some remarkable differences between the data and the 



responsivity at the edges of this structure and also the 
same kind of profile is found in the revl06 data. Therefore, 
we believe it is real. 

HD 45677: The 3 features aroimd 14 /^m, the 13.5, 
13.8 and 14.2 micron features, are clearly present in the 
spectrum. The 20.7 and 21.5 micron feature are rather 
prominent and almost as strong as the 23 micron complex. 

89 Her. The (instrumental) features around 14 /xm are 
all present, but the 20.7, 21.5 and 26.1 micron features are 
not detected. 

MWC 300: The (responsivity) features around 14 jim 
are all present in this source. This is one of the few sources 
where the 21.5 micron feature is not present, while the 
20.7 micron feature is rather prominent. This is evidence 
that these two features are (at least partially) caused by 
different dust species. The 26.1 micron feature is present 
and rather prominent. 

HD 44179: The 13.5 and 14.2 micron features are 
found, but might be due to carbon rich material that also 
peaks at these wavelengths (Hony et al. 2000). There is 
no indication of a 13.8 micron feature. There is a weak 
26.1 micron feature present. In the LWS regime a 48 mi- 
cron feature is found, but it could not be resolved into 
the two separate components. No 91 micron feature could 
be detected because the spectrum was unreliable at these 
wavelengths. 

Roberts 22: Interestingly the 14.2 feature seems present 
but without the 13.6 feature. The 26.1 micron feature 
seems absent and there is also no indication for features 
around 48 /im. The presence of the other feature are uncer- 
tain and will be discussed in the Appendix. There is some 
uncertainty over the presence of a 91 micron feature. 

4.8.2. Outflow Sources 

Vy 2-2: The features around 14 /zm and the 26.1 micron 

feature are not detected. Whether this is due to the noise 
level or that they really are absent is difficult to say. A 
comparison with other stars suggests the flrst possibility. 
While the 20.7 and 21.5 micron features are somewhat 
uncertain, they fit the average profile quite well. Never- 
theless, the noise around these wavelengths makes them 
suspect. 

HD 161796: The 21.5 micron feature is affected by the 
differences seen in the up and down scan. Comparison with 
the rev071 and other data suggests that the 21.5 micron 
feature is probably real, while other features in this wave- 
length region are probably not. The two features around 
48 micron are found blend. 

OH 26.5+0.6: Although all features are in absorption, 
the 3 features around 14 /zm appear in emission. This 
strengthens our suspicion that these features are instru- 
mental. The 20.7 micron feature is thought to be in ab- 
sorption, although it is rather blue-shifted with respect to 
the emission spectra from other sources. The 26.1 micron 
feature might have some self-absorption. 
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HD 179821: There arc indications of a broad feature 
around 13 micron, which is also seen in HD 161796 and 
AFGL 4106. However, it is more likely that it is just con- 
tinuum and that the 18 micron amorphous silicate feature 
already starts around 14.5 fjtm.. The 14 micron (responsiv- 
ity) features were not found in this source. A 20.7 micron 
feature has been found, without a 21.5 micron feature. 
Also the 26.1 micron feature has been found. The 2 fea- 
tures around 48 fjbm. were nicely split and therefore mea- 
sured separately. 

AFGL 4 106: There is a broad feature around 13 /xm 
which is also seen in HD 161796 and HD 179821. It is 
thought to be just continuum and probably induced by an 
early onset of the 18 micron amorphous silicate feature. 
On top of this broad structure, we see the (instrumental) 
features around 14 /im. The 20.7 and 21.5 micron fea- 
tures are found as well as the 26.1 micron feature, which 
is rather broad. We also detect the 48 micron blend. 

NML Cyg: The 13.6 and 14.2 micron features are 
blended with absorption features from HCN at 14.02 /xm 
and an unknown broad absorption at 13.25 /xm. These fea- 
tures are seen in all datasets. The 20.7 and 21.5 micron 
features were detected. The 48 micron blend has also been 
detected (and measured). 

IRC+1042: The features aroimd 14 /im are all found. 
The 20.7 and 21.5 micron feature are also present, al- 
though the first one is the broadest one found in our sam- 
ple. The 26.1 micron feature is weak but present. The 
47.7 and 48.6 micron feature are also found, with the last 
one, in contrast to what has been found for other silicate 
features in this star, the most red-shifted found in our 
sample. 

5. Discussion £ind conclusions 

We have made an extensive study of the infrared spec- 
tra of 0-rich environments of evolved stars. The spectra 
show a remarkably rich variety of features, most of which 
never seen before. We have defined seven complexes at 10, 
18, 23, 28, 33, 40 and 60 /im. Each of these complexes 
is a convolution of several contributions, which vary in- 
dependently, giving each spectrum a unique appearance. 
We have derived average shapes for the complexes. The 
large sample enabled us also to identify weak structures. 
We have quantified the properties of the complexes and of 
the components that make up these complexes in terms 
of peak position, band strength and band width. We find 
that the band strength correlates with the geometry of 
the circumstellar environment, and have accordingly con- 
ducted a systematic separate inventory of the bands and 
complexes for disk sources and outflow sources. We have 
presented a detailed description of the shapes of the bands 
and complexes of the individual stars, and how the spec- 
tra of individual objects differ from the average complex 
shapes. In Paper II, we will discuss, compare and identify 
the average complexes. 



The richness of the spectra undoubtedly contains im- 
portant information related to the conditions in which the 
grains were formed and their thermal history. The crys- 
talline nature of the carriers of most of the narrow bands 
(see also papers II and III) allows a detailed chemical anal- 
ysis of these components, in a way which is unprecedented. 
In this way, it may become possible to trace the evolu- 
tion and processing of grains from their place of birth 
(in the outflows of evolved stars), through the interstel- 
lar medium, into star forming regions and proto-planetary 
disks. For instance, Molster et al. (1999a) suggested that 
in the disk sources low temperature crystallization takes 
place. The crystallization of the silicates in outflow sources 
is likely to have taken place close to the star, i.e. high tem- 
perature crystallization. Therefore the difference found in 
the spectra and thus dust properties might have to do with 
the main crystallization process in these sources. 

The data set presented in this study will serve as a 
starting point for more detailed investigations relating to 
the identification of the carriers of the bands, and for cor- 
relation studies. These will be presented in papers II and 
III. In Paper II we will discuss the average complexes and 
will identify the carriers of the bands using laboratory 
spectra of a variety of plausible components. In Paper III 
we will study the correlations between the wavelengths 
and strengths of the bands, and derive fits to the bands 
using simple radiative transfer and laboratory spectra. 
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6. Appendix; Notes on individual stars 
6.0.3. IRAS09425-6040 

Wc have 2 SWS AOTOl observations for this star, and 
unfortunately no LWS data. We have slightly tilted band 
2A so that the different subbands were nicely connected 
with each other. 

There is C2H2 absorption at 7.8 /um, also HCN ab- 
sorbs around these wavelengths. The change in slope 
around 10.64 /im suggests the presence of a 10.7 mi- 
cron feature. However, no other features are found, which 
makes it somewhat doubtful. The 11.3 SiC feature pre- 
vents the detection of the 11.4 micron feature associated 
with forsterite. But, because of the absence of structure 
at shorter wavelengths, it is likely not present. 

Because of the noise level, the 15.3, 15.9 and 16.2 
micron features are measured as one, but they seem all 
present in this star. There are indications in the spectrum 
that there is a 17.5 micron feature, but this feature has not 
been measured separately. There might be an underlying 
amorphous silicate feature, but it has not been measured. 
The 24.0, 28.8 and 34.1 micron features might be present 
but have not been measured separately. There is a hint 
of the 31.2 micron feature in the rev084 data. However, 
this cannot be confirmed on the basis of the rev254 data 
because of the noise level in this area. 

6.0.4. NGC 6537= HD 312582= Hen 2-340= 
IRAS18021-1950 

For NGC 6537 wc have 2 SWS and 2 LWS spectra. Unfor- 
tunately, the second SWS spectrum was mispointed and 
could hardly be used, even to check uncertain features in 
an independent datasct. 

The origin of the sharp feature at 18.88 /um is un- 
known. A rise in the spectrum is found in both the up 
and the down scan. A similar feature is seen in AC Her, 
but only in one of the two datasets. No other stars show 
this feature, including NGC 6302, which is almost a twin 
to this source. Therefore, we still have some doubts about 
this feature. Independent of the status of this feature, it 
influenced our measurement of the broad 18.9 micron fea- 
ture underneath. The 24.5 micron feature is severely in- 
fluenced by the [Ne V] emission line. The 27.6 and 28.2 
micron features are measured as one. 

6.0.5. NGC 6302= Butterfly Nebula= HD 155520= 
Hen 2-204= IRAS17103-3702= CD-36°11341 

Since NGC 6302 was one of the calibration sources, a lot 
of data was available. We had an AOTOl and a com- 
plete wavelength coverage in AOT06, together with 7 LWS 
datasets. The SWS pointings were slightly offset from the 
centre. Therefore, we multiplied the whole SWS spectrum 
to fit the LWS spectrum. The final spectrum is in agree- 
ment with the IRAS broad band fluxes. 



The 13.6 and 14.2 micron features seem to be present 
in the AOTOl dataset, while in the AOT06 dataset only 
the 14.2 seem to be present. However, strong atomic lines 
at these wavelengths make a clear detection of these fea- 
tures impossible. The spectrum of NGC 6302 shows fea- 
tures at 15.3 and 15.8 /um. However, the strong [Ne HI] line 
in between these features influences them. The difference 
in the up and down scans also made this a problematic 
area. The feature at 16.3 fim feature is only seen in band 
3C of the AOTOl speed 4 dataset, and not in band 3A of 
the same dataset. These last 3 features are also not visi- 
ble in the AOT06 spectrum, therefore we decided not to 
include them. The sharp feature at 19.1 /im is found in 
both the up and down scans. Also the AOT06 spectrum 
shows evidence for a narrow feature at that wavelength. 
The source of this feature is unknown. It seems too narrow 
for a dust feature. We cannot exclude a forbidden emission 
line, although it is quite broad with respect to the other 
lines in this wavelength region and we don't have a clear 
identification yet. The 24.5 micron feature is difficult to 
measure since there is another [Ne V] line close to it. Also 
for the 26.1 micron feature we face the problem of a rather 
strong line on top of it ([O IV]), which prevents us from 
getting a reasonable estimate for this feature. If the 40.5 
micron feature is a blend of 2 features, as is suggested in 
paper II, these two features nicely blend into one structure 
peaking at 40.5 /xm. 

6.0.6. MWC 922= Hen 3-1680= AFGL2132= 

IRAS18184-1302 

The SWS spectrum of MWC 922 is derived from a com- 
bination of 2 AOTOl speed 2 observations. Features were 
treated as real if they were seen in both observations. 

The 9.1 and 9.5 micron absorption features seem to 
form a blend, this might be caused by a known instru- 
mental feature at about 9.35 /im. The 9.8 micron feature 
is very clear and seen in both datasets (up and down scan), 
the same holds for the 10.6 micron feature. The 11.3 mi- 
cron feature, if present, is completely overwhelmed by the 
11.3 micron PAH feature. 

As discussed in the text, the 18 micron complex might 
suffer from interstellar absorption. It was impossible to 
disentangle the circumstellar amorphous silicate feature 
from the interstellar absorption and the circumstellar crys- 
talline silicate emission and therefore it has not been mea- 
sured. The (interstellar) absorption by amorphous silicates 
prevents us from measuring the features around 16 /itm. 
Based on a comparison with other stars, these features 
seem to be present, but there remains uncertainty about 
their strength and FWHM. The 18.0, 18.9 and 19.5 mi- 
cron features are clearly present. The small structure (3 
features) on top of the 18.1 micron feature is located at 
the same position as in HD 44179. The presence of a weak 
feature at 17.43 /im is uncertain, since the two datasets 
do not fully agree around this wavelength. The feature at 
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33.5 /im is the [S III] Hnc. No separate 36.0 micron feature 
can be distinguished in the plateau. The 47.7 and 48.6 mi- 
cron features were difficult to disentangle and therefore 
measured as one. The [O I] line is visible at 63.18 ^m. 
There is a small extension at the long wavelength side 
of the 69.0 micron feature. It might be caused by an un- 
known dust component, or by a population of relatively 
warm forsterite grains, which will emit at a slightly longer 
wavelength (Bowey et al. 2000). 

6.0.7. AC Hcr= HD 170756= SAO 86134= 
IRAS18281+2149 

For AC Her, we have two SWS and one LWS datascts. No 
special treatments were necessary during data reduction. 

The feature at 18.85 micron is only seen in the rev520 
datasct and not in the revl06 data. Although a similar fea- 
ture has been found in NGC 6537, it makes this feature 
questionable. There is a hint in the spectrum of the 22.4 
micron feature and we have measured it, but the reader 
should be aware of the noise level in this part of the spec- 
trum. The strong feature at 22.69 /um looks rather promi- 
nent and seems also to be present in HD 44179, where this 
feature is seen in both the AOTOl as well as the AOT06 
data. However, there is no indication of this feature in 
the rcvl06 spectrum of AC Her. It should be noted that 
the responsivity drops around that wavelength. It is clear 
that this feature should be treated with caution. The re- 
ality of the sharp peak around 31.0 fim is questionable. 
This noise peak sits probably on top of a genuine 31.2 
micron feature, which seems also present in the revl06 
spectrum. The shape of the underlying 31.2 micron fea- 
ture corresponds rather well with the mean shape of the 
31.2 micron feature. The prominent feature at 32.61 /um is 
present in both the up and down scans of rev520 (although 
much more prominent in the up scan). However, nothing 
is seen in the revl06 dataset and we have some doubt 
about its reality. Nevertheless, it should be noted that its 
wavelength position is not atypical for a 32.8 micron fea- 
ture. The 33.7 and 34.1 micron features are measured to- 
gether. In the rcv520 datasct, there is an indication of the 
presence of a 35.9 micron feature and the absence of the 
34.9 micron feature. However, it seems vice versa for the 
rcvl06 data, where a 34.9 micron feature seems present 
but hardly anything at 35.9 /um. It is clear that the noise 
level in both datasets play tricks on us here. We interpret 
this as evidence for the presence of the 35 micron plateau 
although it is not clear where this plateau ends. 

6.0.8. HD 45677= FS CMa= MWC 142= SAO 151534= 
AFGL 5195= IRAS06259-1301 

For HD 45677 we only had one SWS dataset and unfor- 
tunately no LWS data. Band 2 B has been rotated to let 
all the bands nicely match. 



In general, the features of HD 45677 are rather sharp 
compared to the other stars. The 34.1 and 34.9 micron 
features are both present but blend in this object and are 
therefore measured as one. The same applies for the 35.9 
and 36.5 micron features. The sharp line at 34.81 /xm is 
due to [Si II]. 

6.0.9. 89 Her= V441 Her= BD-h26°3120= HD 163506= 
IRAS17534-t-2603= AFGL 2028= SAO 85545 

For 89 Her, we also have redundancy in the SWS domain 
with two datasets. For the LWS wavelength range, we have 
only one dataset. To match the different bands with each 
other, we are forced to slightly tilt band 2B. 

There is uncertainty about the 15.2 micron feature. 
The structure seen at these wavelengths seems real, be- 
cause it is found in both the rev082 and the rev518 data. 
However, it is rather blue-shifted and also very broad. 
Therefore, the question arises whether the feature found 
at 15.0 /xm in this object is the same as the 15.2 micron 
feature. 

We have fitted the smooth 23 micron complex with 
3 features but these were not well constrained. Therefore 
caution should be taken in interpreting the results of this 
complex, even though the errors might look small. 

The 31.2 micron feature is not detected. The structure 
that is seen around this wavelength is noise. 

6.0.10. MWC 300= V431 Sct= AFGL 2170= 
IRAS18267-0606 

For MWC 300 we unfortunately have only one SWS spec- 
trum and no long wavelength data. The substructure 
found in the 10 micron complex, which might be at- 
tributed to crystalline silicates, is always seen in only one 
scan and we therefore do not trust its appearance. We have 
therefore not measured them. There might be structure 
around 16 /zm, but this is severely blended by interstellar 
silicate absorption and has therefore not been measured. 
The amorphous silicate feature detected at 17.03 /zm is 
probably due to interstellar absorption and not a cir- 
cumstellar amorphous silicate emission feature, which has 
therefore not been measured. The up and down scans dif- 
fer for the sharp peak at 28.68 /zm; in one it is present, 
but in the other it is not. But the drop aroimd 28.9 /im is 
seen in both the up and down scan, confirming the pres- 
ence of the 28.8 micron feature. There is an indication of 
the 32.2 micron feature but at a very low level. Therefore 
it has not been measured. The up and down scans do not 
correspond at 42.5 micron. Therefore, we have removed 
this feature before measuring the 40 micron complex. The 
shallow slope of the blue side of the 40.5 micron feature 
suggests the presence of a 39.8 micron feature. The pres- 
ence of the 44.7 micron feature can neither be proven nor 
rejected. 
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6.0.11. HD 44179= Red Rcctangle= BD-10°1476=: 
IRAS06176-1036= AFGL 915= SAO 151362 

For this source wc have an AOTOl spectrum and some 
parts done with a higher resolution (AOT06). Also an 
LWS spectrum is available. 

The feature at 17.07 /xm is likely to be a blend of the 
16.7 and 17.0 micron features. The reality of the 20.7 
micron feature is uncertain since it is very broad in the 
AOTOl spectrum (although seen in both the up and down 
scans), but nothing is seen in the AOT06 spectrum. There- 
fore, caution should be taken in interpreting this feature 
in this star. The structure around 23.2 /im is not found 
in all datasets and differs between the up and the down 
scans. Therefore, wc masked this part before measuring 
the complex. Comparing the 28 micron complex with 89 
Her suggests that the 26.8 micron feature is also present 
on the blue side of the 27.6 micron feature. 

6.0.12. Roberts 22= Hen 3-404= Wray 15-549= 

AFGL 4104= OH284. 18-0.79= IRAS10197-5750 

When ISO observed Roberts 22 with the SWS it was off- 
set, and a lot of flux was missed. Increasing apertures gave 
also flux jumps between the different subbands. In total, 
we have 2 SWS and 1 LWS spectrum for this source. 

We have measured the 18 micron amorphous silicate 
band. However, it is likely that the 18.0, 18.9 and 19.5 mi- 
cron features are present and are therefore included in this 
measurement, but the noise level prevents firm conclusions 
about their presence. The 15.9 and 16.2 micron features 
are measured together. The 23 micron complex has been 
measured as one feature, since it was too noisy to separate 
the different features. The 28 micron complex suffers from 
flux jumps between the different bands. Therefore, there is 
some uncertainty about the continuum in band 3E. Still, 
we were able to measure the features in this band since 
we used a local continuum for these measurements. How- 
ever their relative strength is questionable. All 3 features 
found in band 3E were seen in both the rev084 as well as 
the rev254 datasets. Before measuring the 33 micron com- 
plex, we masked the narrow (0.1 /zm) peak at 33.3 /zm, 
since it was only found in one scan direction and not in 
the rev084 data. The 35.9 -|- 36.5 micron features were 
measured as one, alike to the 39.8 + 40.5 micron features. 
The 20.7 and 21.5 micron features are found in rcv254, 
where their shape resembles the average shape of these 
features quite well. However, because nothing is found in 
the rev084 dataset, and because of the low S/N in this 
part of the spectrum, these features are suspect. There is 
some concern about the data beyond 80 /xm. Therefore, 
the detection of the 90 micron feature in Roberts 22 is 
uncertain. 



6.0.13. Vy 2-2= IRAS 192 19+0947 

Vy 2-2 is the other source on which ISO was not centered 
correctly. This resulted in flux jumps and intensity losses. 
We had one SWS and one LWS spectrum. The whole SWS 
spectrum suffered from the mispointing of the satellite. 
This reduced the flux level and therefore the signal to 
noise ratio. 

The reality of the feature at 11.6 /im is severely 
doubted, since the up and down scans do not agree with 
each other. It is therefore a good indicator of the noise- 
level in that part of the spectrum. The amorphous silicate 
features, clearly present in the 10 micron complex, could 
not be detected at the 18 micron complex, likely due to 
the noise level. 

The feature at 25.1 fim is too broad to be a forbidden 
emission line. However it also seems narrower than the 
normal 25.0 micron feature. The noise in this part of the 
spectrum is quite severe, therefore we have some doubts 
about this feature. If other data confirm this to be the 

25.0 micron feature, it would be the strongest 25.0 micron 
feature in the sample. The overall intensity from 24.2 to 
25.2 /xm suggests that the plateau is present and this has 
been measured (excluding the uncertain strong peak at 

25.1 Atm). 

There were a lot of spurious features in band 4, which 
were removed by hand during the data reduction. There 
are hints in both the up and down scans for the different 
features in the 35 micron plateau, therefore we have mea- 
sured them. However, the signal to noise ratio in this re- 
gion makes these features somewhat uncertain and there- 
fore no firm conclusions should be drawn based on these 
features alone. Also the 40 micron complex suffers from 
the high noise level. The features seen in the spectrum 
provided by the ISO-data archive in Vilspa are probably 
due to glitches and glitch tails, which were found in this 
part of the spectrum. We have removed these spurious 
points by hand. We have tried to split the 40 micron com- 
plex into a 40.5 micron feature and a 43.0 + 43.8 micron 
blend and measured the result. 

6.0.14. HD 161796 = V814 Her= BD+50°2457= 

IRAS17436+5003= AFGL 5384= SAO 30548 

We got a nice dataset for HD 161796, including an AOT06 
observation of band 4. In the final spectrum we attached 
this AOT06 spectrum between 29.03 and 43.17 fim. 

A broad feature around 13 micron may be present, 
which is also seen in AFGL 4106 and HD 179821, however 
in the latter spectra it resembles continuum with an early 
rise of the 18 micron silicate feature, starting at about 
14.5 /im. The 18 micron complex is located on the slope 
of a steep spectrum and there is some uncertainty about 
the strength of the 18 micron feature. The up and down 
scans do not completely agree around this feature. Still 
the peak position and FWHM agree very well with the 
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feature seen in the mean outflow spectrum. So wc conchidc 
that it is likely that there is a feature at this wavelength, 
but its strength is somewhat controversial. Because of the 
presence of the crystalline silicate features, it is more diffi- 
cult to measure the amorphous silicates. This dust feature 
also hampers from our method of measurement. We fitted 
Gaussians while this feature is not a Gaussian. The small 
feature which is seen at 17.5 /xm is noise. The shape of the 
20.5 micron feature is similar to the shape in the mean 
spectrum, however its strength is similar to the features 
at 22.0 and 22.4 /xm, for which the up and down scans 
differ significantly. Therefore, the detection of the 20.5 mi- 
cron feature is only marginal. Because of the reasonable 
uncertainty about the 22.4 micron feature, it has not been 
measured. The 28.8 micron feature might be present but 
one should take into account that we only measured it af- 
ter masking the spurious feature at 29.2 /xm. It should be 
noted that the 30.6 micron feature is found in the AOT06 
data and in the rev342 data but nothing is seen in the 
revOTl dataset. The spectrum between 29.03 and 43.17 
jum is from the AOT06 observation. The 35.9 and 36.5 
micron features have been measured as a blend. The 35.1 
micron feature is probably somewhat underestimated and 
blue-shifted. The status of the feature at 44.09 /xm is not 
completely sorted out. It is seen in all data sets, but the 
up and down scans do often differ at this point. For now, 
we assume that it is the result of a blend of the 43.8 and 
the 44.7 micron features. 

6.0.15. OH 26.5-^0.6= V347 Sct= IRAS18348-0526= 
AFGL 2205 

For OH 26.5+0.6 we got 2 LWS and one SWS dataset. 

Up to 28 /im all features are in absorption. The amor- 
phous 10 micron feature suffers from a high degree of ab- 
sorption and no Gaussian profile is fitted to it, the centre 
of the absorption is at 9.9 ± 0.1 /im. The substructure in 
the 10 micron absorption band might be real, it is seen 
in both the up and down scans, but the positions do not 
match with other sources in our sample. The 23.0 and 23.7 
micron features are measured as a blend. It is not com- 
pletely clear whether the 27.6 and 28.2 micron features are 
present or not. They are located at the border where the 
spectrum goes from absorption to emission. It seems that 
there is some excess emission at these wavelengths, but 
no clear peak has been seen. This suggest that they suffer 
from self absorption, which makes it impossible to measure 
them with Gaussians. So we have not measured these fea- 
tures, although they seem to be present. The 33 micron 
complex without the plateau has been measured as one 
feature, since it was not possible to make sure how to sep- 
arate the different features, for which are indications for 
their presence. Although there is an indication of 2 differ- 
ent features in the 60 micron complex of OH 26.5-1-0.6 we 
have measured it as one. There seems to be self- absorption 
in the 26.1 micron feature, which hampers the measure- 



ments of this feature. The 48 micron blend has again been 
measured as one due to the noise level in that part of the 
spectrum. 

6.0.16. HD 179821 = V1427 Aql= BD-00°3679= 

IRAS191 14-^0002= AFGL 2343= SAO 124414 

Our data products of HD 179821 consist of two SWS and 
one LWS spectrum. 

The 18 micron complex lies on a very steep slope, 
which makes it difficult to identify and measure the dif- 
ferent features. The 18 micron amorphous silicate feature 
might be about 10% too strong since it has been mea- 
sured including the crystalline silicate features. The 15.9 
and 16.2 micron features were measured together, since it 
appeared very difficult to separate them. The 17.0 micron 
feature is very weak, but its width corresponds well with 
the mean feature shape. The detection of the 18.9 micron 
feature is uncertain, it has therefore been measured to- 
gether with the 18.0 micron feature. The 24.5 and 25.1 
micron features have been measured as one. The 40.5 mi- 
cron feature shows a relatively shallow slope towards the 
shorter wavelengths, this is fitted with an extra Gaussian, 
which is attributed to the 39.8 micron features. This ex- 
tra Gaussian influences the strength and probably also the 
derived width of the 40.5 micron feature. The sharp drop 
around 67 /im might have to do with the border between 
two different bands in the LWS spectrum. We remind the 
reader that we have not performed any shift of the sub- 
bands for this LWS spectrum. 

6.0.17. AFGL 4106= HD 302821= CPD-58°2154= 
IRAS10215-5916 

There are three SWS datasets for AFGL 4106 and one 
LWS dataset. We had to rotate band 3E to match band 
3D and 4. 

The 18 micron amorphous silicate feature might be 
about 10% - 20% too strong, since it has been measured 
including the crystalline silicate features. Again some fea- 
tures have been measured as a blend, such as the 15.9 -I- 
16.2 micron features, the 16.7 -|- 17.0 micron features and 
the 47.7 -I- 48.6 micron blend. There are slight indications 
of a 22.4 micron feature, however we decided not to mea- 
sure it because of the noise level, which makes this feature 
rather uncertain. 

6.0.18. NML Cyg= V1489 Cyg= IRC-|-40448= 
OH80.8-1.9= AFGL 2650 

For NML Cyg there was a nice dataset available, with 4 
SWS, of which one was an AOT06 observation of band 
4, and 3 LWS observations. We have included this high 
resolution spectrum in our flnal spectrum. Band 2A has 
been tilted to fit the other bands. 
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The 8.3 micron feature might be present, but is diffi- 
cult to separate from the amorphous silicate feature which 
still peaks at these wavelengths. The structure around 
10.07 fim coincides with a feature in the responsivity and 
is therefore not trusted. The structures around 11.05 /um 
are due to the apparently not completely correctly re- 
moved 11.1 micron responsivity feature. 

The 13.6 and 14.2 micron features are influenced by 
molecular absorption features at 13.25 and 14.02 /zm. 

We assumed that there was no 13.8 micron feature and 
that therefore the flux at that wavelength was continuum. 
The features were measured with respect to that contin- 
uum. 

The 16.9 and 17.2 micron features are measured as 
one. If there is a very weak 19.5 micron feature it has 
been measured together with the 18.9 micron feature. At 
14.96 /im gas phase CO2 absorption has been detected. 

The 18 micron amorphous silicate could not be ex- 
tracted from the data, it is probably between emission and 
self absorption. The 24.5 micron feature is rather blue- 
shifted and strong compared to the other 24.5 micron fea- 
tures. It gives NML Cyg the broadest appearance for this 
complex of all sources. We cannot exclude that this feature 
is not the same as the other 24.5 micron features. 

The end of band 3D and the beginning of Band 3E are 
not OK due to leakage, therefore no features are measured 
in this wavelength range There might be a 39.8 micron 
feature present, however this is difficult to measure with 
so many water lines present in that region. Therefore we 
measured it as one feature, neglecting the presence of the 
water lines. This might lead to an overestimation of the 
flux level of the features. 

The 47.7 and 48.6 micron features were measured as a 
blend. 

6.0.19. IRC+10420= V1302 Aql= IRAS19244-|-1115= 
AFGL 2390 

The last source in our sample is IRC+10420. We got 
complete spectral coverage in the SWS and LWS wave- 
length range. Part of the SWS-spectrum was done with 
an AOT06 which is inchided in our final spectrum. The 
infrared spectrum of IRC+10420 shows some structure in 
the 10 micron complex which lines up quite nicely with 
crystalline silicates in other sources. The weak 10.6 mi- 
cron feature is uncertain, since it is not consistent with 
the feature seen in the mean disk spectrum feature. The 
11.3 micron feature is consistent with a feature seen in AC 
Her but very weak. 

The 16.2 micron feature is somewhat unsure, since it 
is clearly seen in band 3c (up and down scan), but not 
in band 3 A (neither in the up nor in the down scan). 
The amorphous feature has been measured without the 
crystalline features and the strength is therefore somewhat 
overestimated, however this effect is negligible compared 



to the uncertainty in contimmm level. The 18.0 -|- 18.9 
micron features are measured as one. 

Band 3D was not corrected for leakage problems in the 
calibration sources, because the flux level in IRC+10420 
is so high that it is very likely that also here leakage has 
its influence. Band 3E is tilted to flt the AOT06 and LWS 
on one side and the rest of the data at the other side. 
The 20.7 micron feature is the broadest one found in our 
sample, still the data seem valid, since the up and down 
scans show exactly the same trend. 

The up and down scans of the AOT06 spectrum differ 
severely between 29.3 and 30 /zm, probably due to memory 
effects for this high flux source. Therefore we have removed 
this part of the spectrum, since also the AOTOl suffered 
from the same problems in this part of the spectrum. 

There is a possibility of a 39.8 micron feature however 
the different datasets contradict each other around these 
wavelengths. Therefore we have measured this feature as 
a 40.5 micron feature. 



